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1.0  INTRODUCTION 


1.1  PROGRAM  GOALS 

The  aim  of  this  program  was  to  study  the  molecular  and 
electronic  structure  of  Trinitrotoluene  (TNT)  in  its  various 
crystal  phases.  Sample  of  cast  TNT,  which  normally  contain  a 
mixture  of  orthorhombic  and  monoclinic  phases,  will  often  develop 
cracks  and  voids  when  orthorhombic  monocrystals  transform  to  the 

denser  monoclinic  phase  upon  annealing.  The  desire  to  learn  more 
about  the  molecular  structure  and  molecular  dynamics  of  TNT  in 
these  two  phases  provides  part  of  the  rational  for  this  study. 
Nitrogen-14  nuclear  quadrupole  resonance  (NQR)  techniques  were 
proposed  as  singularly  well  suited  to  the  problem  and  because  they 
complement  x-ray  and  other  crystallographic  studies.  It  was 
further  proposed  to  study  the  molecular  dynamics  of  TNT  in  the 
solid  state.  Finally,  explosive  systems,  such  as  RDX  and  HMX, 
were  proposed  as  candidates  for  an  NQR  study. 

1.2  SUMMARY  OF  PROGRAM  RESULTS 

We  have  established  that  NQR  studies  are  well  suited 
to  the  characterization  of  explosive  substances  such  as 
TNT,  RDX,  and  HMX. 
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The  complete  Nitrogen-14  NQR  Spectrum  of  TNT  was  obtained 
for  both  the  stable  phase  (monoclinic)  and  the  unstable 
(orthorhombic)  phase.  It  was  established  that  quadrupole  coupling 
data  are  very  sensitive  to  and  can  yield  a  measure  of  the  angle 
of  twist,  8,  of  the  plane  of  the  N02  group  away  from  the  plane  of 
the  benzene  ring  about  the  CN  bond.  The  temperature  dependence 
of  the  NQR  data  was  then  used  to  obtain  <A0>  RMS,  the  root  mean 
square  angular  displacements  about  the  equilibrium  value  of  Q 
at  room  temperature.  Large  values  of  <A0>  RMS  were  found  in  excess 
of  10°  at  room  temperature.  Also,  large  variations  in  <A0>  RMS 
were  observed  among  different  N02  groups. 

For  the  high  explosive  1 , 3 , 5 , 7- tetramtro-1 , 3 , 5 , 7-tetraaza- 
cyclooctane  (HMX) ,  an  analysis  of  the  Nitrogen-14  NQR  spectrum 
and  its  temperature  dependence  yield  details  of  the  motion  of  the 
molecule  in  the  solid.  The  Nitrogen-14  NQR  spectrum  of  RDX  at 
77°K  was  also  reported. 

Finally,  two  instrumentation  advances  were  made.  A  new  pulse 
sequence,  labelled  SORC  for  Strong-Of f-Resonant  Comb  of  RF  pulses, 
was  introduced  and  some  of  its  properties  were  studied.  In  addition, 
moved  by  a  desire  to  display  as  much  as  possible  of  the  crowded 
NQR  spectrum  of  TNT  as  possible,  Fourier  transform  (FT)  techniques 
were  extended  to  obtain  close  to  an  order-of-magnituae  increase 
in  the  breadth  of  the  spectrum  displayed  in  a  single  Fourier 
transform. 
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2 . 0  EXPERIMENTAL  TECHNIQUES 
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Of  the  many  experimental  methods  suitable  for  the 
detection  of  nuclear  quadrupole  coupling  in  solids,  pulsed  NQR  in 

4 

conjunction  with  Fourier  transform  techniques  offer  important 
advantages  for  the  current  project.  Principal  among  these  is  the 
high  resolution  with  which  individual  line  shapes  can  be  detected 
and  displayed.  This  enables  us  to  resolve,  for  example,  all  the 
features  of  the  remarkably  crowded  Nitrogen-14  NQR  spectrum  of 
a-TNT:  twenty-four  lines  were  detected  in  the  frequency  region 
700-900  kHz. 


A  brief  discussion  of  the  basic  concepts  of  the  theory  of  NQR 

is  presented  in  Appendix  A.  The  pulsed  NQR  spectrometer  used  in 

this  work  was  developed  at  Block  Engineering  and  has  been 
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published  elsewhere  .  In  the  next  sections  we  will  describe  some 
of  the  specific  advances  in  experimental  techniques  that  were  made 
during  this  project. 

2.1  PULSED  TECHNIQUES  FOR  NITROGEN -14  NQR 


Extensive  signal  enhancement  is  often  required  before 
Nitrogen-14  nuclear  quadrupcle  resonance  lines  can  be  displayed 
with  satisfactory  signal-to-noise  ratio.  The  pulse  sequence  most 
suited  to  a  particular  situation  and  the  data  handling  used  are 
chosen  with  this  in  mind.  Four  possible  pulse  sequences  are 
discussed  below.  The  optimum  choice  depends  on  the  value  of  the 
relaxation  times  of  the  particular  transition  under  study. 
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2.1.1  Free  Induction  Deca' 


The  NQR  response  of  a  single  crystal  to  a  resonant  pulse  of 
RF  irradiation  is  completely  analogous  to  the  spin  1/2  NMR  case. 
One  difference  is  that  the  intensity  of  the  free  induction  decay 
(FID)  response  depends  on  the  orientation  of  the  radio-frequency 
field  with  respect  to  the  electric  field  gradient  (EFG) 
principal  axes  system.  For  a  spin  1=1  system,  the  three 
resonance  lines  v  +  ,  v_  and  vd  can  be  observed  only  when  the 
irradiation  field  is  oriented,  respectively,  along  the  x,y,  and 
z  principal  axes  of  the  EFG  tensor.  As  an  example,  for  a  v_  line, 
the  expected  value  of  the  lacsnetization  along  the  y  axis  is 
proportional  to: 


sin  (/2  ant  J  cos  (2irv  t) 

1  w  - 


(2. 1.1-1) 


where  ^  =  Y measures  the  intensity  of  the  irradiation,  and  ty 
is  the  duration  of  the  irradiating  pulse.  The  maximum  response  is 
obtained  when: 


/2  on  t  =  /2  Y  H.t  =  *  /  2 

I  w  r  w 

in  analogy  with  the  NMR  case. 


(2.1. 1-2) 


In  NQR,  however,  the  sample  often  consists  of  a  polycrystalline 

powder.  In  such  cases  a  convolution  for  all  orientations 

must  be  made.  The  result®  is  that  the  sin  (/2w1t  )  function  in 

the  expression  for  the  expectation  value  of  the  magnetization 

becomes  a  Bessel  function,  J^(/2u>^tw).  This  function  has  its 

first  maximum,  analogous  to  a  "90°  pulse",  at  /2w.t  =  0.66tt  and 

i  w 

not  at  0.5 *  as  the  sine  function.  Similarly,  the  f-rst  null, 
corresponding  to  a  "180°  pulse"  occurs  for  a  value  of  the  argument 
equal  to  1.43*  rather  than  simply  to  a  value  of  n. 
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The  time  constant  of  the  free  induction  decay  is  labelled 

T2*.  In  solids  this  parameter  is  usually  small,  given  the 

rigid-lattice  broadening  of  energy  levels  caused  by  the  magnetic 

7 

dipole-diploe  interaction.  Therefore,  FID  signals  are  usually  at 

least  partially  masked  by  the  instrument  dead-time  or  recovery-time 

RF 

following  each  intense  pulse.  For  these  reasons  FID  signals 

14 

are  often  not  observed.  This  is  so  even  in  the  case  of  N  NQR 
where  the  magnetic  dipole-dipole  interaction  is  quenched  to  first 

g 

order  because  the  nuclear  spins  have  integer  values.  . 

2.1.2  Spin  Echo  and  Carr-Purcell  (CP)  Sequence 

Following  a  FID  experiment  in  a  case  where  T2  (spin-spin 

relaxation  time)  >>T2*  (FID  time  constant)  it  is  possible  to 

recall  part  of  the  magnetization  not  lost  through  T2  processes  by 

applying  a  "ISO0  pulse"  at  a  time  x  after  the  first  pulse.  As  is 

well  known,  an  "echo"  will  form  at  the  time  2x;  this  echo  signal 

9 

can  be  repeatedly  recalled  at  integral  multiples  of  this  time, 

2nx,  by  the  application  of  additional  "180°  pulses”  at  times 
(2n-i)x.  The  amplitude  of  the  resultant  echo  train  decays  with 
time  constant  T2  (Figure  2. 1.1-1).  Signal-to-noise  enhancement 
can  be  obtained  by  coherently  adding  successive  echoes  in  the 
sequence.  The  optimum  time  for  co-addition  is  easily  shown 
to  be  1.26t2.  Due  to  the  generally  small  value  of  T2  in  solids, 
however,  this  method  does  not  result  in  an  appreciable  enhancement 
of  the  signal-to-noise  ratio. 
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Figure  2.1 .1-1.  Parameters  for  Carr-Purcell  Sequence 


A  major  advance  in  S/N  enhancement  was  made  by  Marino  and 

11  12 
Klainer  when  an  adaptation  of  the  Ostrof f-Waugh  sequence  was 

applied  to  NQR.  The  sequence  is  essentially  a  Meiboom-Gill- 

modified-CP.  sequence^-0  where  all  the  pulses  have  the 

same  flip  angle  of  "90°",  which  means  0.66 t  in  the  NQR  of  an  I  =  1 

nucleus.  In  Figure  2. 1.3-1  (NaN02  at  77°K,  at  resonance)  and  Figure 

2. 1.3-2  (NaNO^  at  77 °K,  slightly  off  resonance)  it  can  be  seen  that 

the  spin  echo  train  in  this  sequence  persists  for  times  of  an  order  of 

the  spin  lattice  relaxation  time,  T^,  and  not  the  much  shorter  T2« 

Coherent  addition  of  the  echoes  in  this  case  results  in 

considerable  enhancement  of  the  S/N,  since  >>  T2  is  the  typical 

situation  in  these  solids.  Marino  and  Klainer  showed  that  the 

1/2 

optimum  enhancement  is  0.64  (T2  /2 i)  where  T2  is  the  effective 

decay  constant  of  the  echo  train  and  2x  is  the  spacing  between 

echoes,  or  equivalently  the  spacing  between  pulses  of  the 

excitation  sequence.  It  was  further  shown  that  the  NQR  effect  was 

13 

completely  analogous  to  the  spin  1/2 ,  NMR  case  discussed  by  Waugh 

in  that  the  decay  constant  T2£  tends  to  as  t  is  reduced  to 

values  less  than  T2«  Furthermore,  for  intermediate  values  of  T2  , 

this  parameter  is  proportional  to  x  ^ ,  again  in  analogy  to  NMR. 

Figure  2.1. 3-3  shows  this  functional  dependence  for  NaNO-  at  77° K. 

14  ^ 

Recently,  Cantor  and  Waugh  have  developed  a  theory  to  explain 
the  main  features  of  this  NQR  effect  using  a  model  of  a 
polycrystalline  solid  with  each  nitrogen  site  having  one  nearest 
neighbor. 

2.1.4  Strong  Off-Resonance  Comb  (SORC)  of  RF  Pulses 

Recently  we  developed  a  new  pulsed  NQR  experiment  which  can 
have  considerable  advantages  in  enhancing  the  S/N  ratio  of  weak 

4 

lines. 
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Fiqure  2. 1.3-1.  Spin-Locked-Spin-Echo  (SLSE)  train  for  v_  line  of  NaN02  at 

77 °K.  Transmitter  frequency  set  at  center  of  resonance  line. 


T'"e  steady  state  response  of  an  ensemble  of  nuclear  spins 
(1=1/2)  in  high  magnetic  field  HQ,  to  a  strong  radio-frequency 
field  H1,  applied  off  resonance  by  Af,  has  long  been  known.8  When 
all  the  conditions  for  the  establishment  of  a  soin  temperature  in 
the  rotating  frame  are  met,  the  x-component  of  the  magnetization,  Mx, 
which  is  experimentally  observable,  is  given  by  the  expression: 


M, 


E1  (2TTAf/y) 

Hl2  +  Hloc  +  <2lTAf/Y>2 


(2. 1.4-1) 


where  Mq  is  the  equilibrium  longitudinal  magnetization,  Y  is  the 
magnetogyric  ratio  of  the  nucleus,  and  H1qc  is  a  measure  of  the 
local  field  at  the  nuclear  site  due  tc  its  neighbors.  Results 
analogous  to  Equation  2.1. 4-1  have  also  been  derived  and  observed 
for  a  quadrupolar  system15  with  nuclear  spin  I  =  3/2  when 
subjected  to  the  same  strong,  long,  off-resonant  irradiation,  H1 . 

The  preliminary  results  obtained,  when  the  irradiation  field 
H1  is  applied  in  a  long  train  of  equally-spaced  identical  pulses, 
are  presented  here.  Although  the  SORC  experimental  data  reported 
here  are  for  a  quadrupolar,  T  =  1  system,  analogous  effects  in  a 
magnetic  system  or  a  quadrupolar  system  with  spin  different  from 
unity  can  be  expected. 

Figure  2.1.4-la  defines  the  parameters  of  the  SORC  sequence. 
Here  a  train  of  radio-frequency  pulses  of  duration  t  ,  spacing  x, 
A f  away  from  exact  resonance  is  applied  to  a  pure  nuclear 
electric  quadrupole  system  in  zero  external  magnetic  field. 


BLOCK 


2-9 


Timing  diagram  for  SORC  sequence  {Strong 
Off-Resonance  Comb  of  RF  pulses) . 


Overlay  of  nuclear  induction  signals 
obtained  in  SORC  observation  window  when 
Af  =  3/2  kHz  and  Af  =  3.5/2  kHz.  Note 
destructive  interference  in  former  case. 
All  data  were  taken  on  the  v  line  of 
NaN02  at  77°K. 

Plot  of  signal  height  obtained  in  the 
middle  and  right-hand  side  of  observation 
window  using  the  SORC  sequence. 
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The  variation  of  the  signal  amplitude  vs.  Af,  the  distance 
from  exact  resonance,  shows  two  features  as  depicted  in  Figures 
2.1.4-lb  and  2.1.4-lc.  First,  the  signal  amplitude  is  modulated  by 
sinusoid  of  period  1/t ,  the  reciprocal  of  the  pulse  repetition 
rate.  This  phenomenon  is  best  understood  by  considering  that  the 
Fourier  transform  of  the  transmitter  pulses  has  periodicity  1/t. 
This  leads  to  successive  maxima  and  minima  in  the  NQR  signal 
(Figures  2.1.4-lb  and  2.1.4-lc)  as  the  transmitter  frequency  is 
changed,  i.e.,  Af  is  varied.  Alternatively,  and  more  naively, 
this  modulation  can  be  interpreted  as  the  destructive  interfernece 
type  I  signals  (FID-like  signals  immediately  following  the  RF 
oulse)  and  type  II  signals  (immediately  preceeding  the  RF  pulse) 
in  theii  overlap  region  as  the  frequency,  Af,  is  slowly  varied. 

Another  featux^  shown  in  Figure  2.1.4-lc  is  the  shape  of  the 

envelope,  possibly  conforming  to  a  function  of  the  type 

2  2 

Af/(A  +  Af  ) ,  such  as  Equation  2 . 1 . 1-2 .  Insufficient  data  have 
been  taken  so  far  to  ascertain  the  degree  of  agreement  with  theory 
on  this  last  point. 

The  nuclear  induction  signals  present  in  the  observation 
window  between  successive  pulses  of  the  SORC  sequence  are  shown  in 
Figure  2. 1.4-2  as  a  function  of  the  pulse  separation t  .  All  data 
were  taken  on  the  v_  line  of  NaNO^  at  77 °K.  The  magnitude  of  the 
type  I  signal  is  then  plotted  vs.  t  in  Figure  2. 1,4-3.  Note  that 
for  msec,  the  magnitude  of  type  I  signals  increases  with  t  ,  as 

might  be  expected  for  a  FID  signal  subject  to  spin-lattice 
relaxation.  On  the  other  hand,  for  ?s.5  msec,  signals  at  both  ends 
of  the  observation  window  are  of  a  comparable  size  and  they  grow 
exponentially  with  decreasing  t  ;  this  is  the  region  of  interest. 
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RELATIVE  SIGNAL  INTENSITY  (ARBITRARY  UNITS) 


t  (MSEC) 


Figure  2. 1.4-3.  A  semi logarithmic  plot  of  the  type  I 

signal  height  (see  text)  vs.  the  pulse 
spacinq  parameter,".  The  data  are 
taken  from  Figure  2.1-6. 
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Figure  2. 1.4-4  shows  the  variation  of  the  type  II  signal  vs. 

<  > for  t  =  3  msec.  The  experimental  points  are  obtained  for 

four  different  values  of  the  instantaneous  field^H^,  obtained  by 
changing  the  pulse  width,  t  ,  at  constant  "flip  angle",  /2  y  -WH1 
The  dotted  line  is  the  curve  F  =  <H^>/  [<H^>2  +  ]  with  B  =  0.05 

G.  The  fact  that  <H^>,  the  average  value,  rather  than  ,  the 
peak  value,  is  the  important  parameter  and  that  there  is  good 
agreement  of  the  data  with  the  form  of  Equation  2. 1.1-2  is  strong 
evidence  that  the  ensemble  of  spins  is  responding  to  the 
time-averaqe  field  of  the  SORC  sequence  in  a  manner  analogous  to 
the  conventional  long,  strong,  off-resonant  pulse. 

The  size  of  the  parameter  B  is  found  to  be  about  0.05G  which 
is  approximately  two  orders  of  magnitude  too  small  for  the  value 
expected  from  the  contribution  of  Af  to  Equation  2. 1.1-2.  This 
discrepancy  is  reduced  by  a  factor  of  5  when  the  experiment  is 
repeated  for  x  =  1  msec,  as  shown  in  Figure  2. 1.4-5.  Comparison  of 

Figures  2. 1.4-4  and  2. 1.4-5  suggests  that  the  pulsed  nature  of  the 
experiment  is  still  very  important  for  t  =  1  msec  and  that 
quantitative  agreement  cannot  be  expected  until  x  is  reduced 
further. 

The  potential  of  this  technique  for  sensitive  detection  of 
NQR  signals  appears  great  since  signals  can  be  obtained  and 
signal  averaged,  at  essentially  100%  duty  cycle.  It  also  appears 
experimentally  that  the  conditions  for  efficient  signal  collection 
on  t/T2,  the  ratio  of  the  pulse  spacing  to  the  spin-spin 
relaxation  time,  are  slightly  less  stringent  here  then  in  the  case 
of  SI,SE.  However,  further  experimentation  is  needed  to  completely 
understand  the  operational  parameters  of  SORC. 
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SIGNAL  HEIGHT  (ARBITRARY  UNITS) 


0.2 
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Figure  2. 1.4-4. 


<HX>G 


A  plot  of  magnitude  of  type  II  signals 
(see  text)  vs,  the  time  average  value 
if  the  KF  field,  <H^>,  using  a  SORC 
sequence.  All  data  were  taken  on  the 
v_  line  of  NaN02  at  77°K  with  = 
3.5/2  kHz;  t  =  3  msec. 
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(ARBITRARY  UNIT) 


Figure  2. 1.4-5. 


A  plot  of  the  magnitude  of  Type  II 
signals  (see  text)  vs.  the  time- 
average  value  of  the  RF  field, 

<H^> ,  using  a  SORC  sequence.  Data 
were  taken  on  the  line  of  NaN02 
at  77°K  with  f  =>  1.5/2  kHz; 
t  =  1  msec. 
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2.2  FOURIER  TRANSFORM  TECHNIQUES  FOR  NITROGEN-14  NQR 


1 


There  are  three  main  reasons  for  doing  FT  spectroscopy: 

a.  Enhancement  of  signal-to-noise  over  CW  methods.  The 
enhancement  is  given  by  the  square  root  of  the  ratio 
of  the  total  width  of  the  spectrum  to  the  typical  line 
width. 

b.  Pulsed  methods  are  singularly  well-suited  to  data 
processing. 

c.  The  line  shape  is  readily  obtained  directly  from  the 
output . 

Of  these,  the  first  reason,  probably  the  most  important  in  other 
disciplines,  is  not  nearly  as  important  in  NQR,  because  the  ratio 
defined  above  is  not  large  and  can  often  be  close  to  unity.  This 
is  so  because  in  solids,  line  widths  are  relatively  large  (a  few 
kHz)  while  the  bandwidths  which  can  be  suitably  irradiated  are  in 
the  10-100  kHz  range.  This  implies  that  the  NQR  spectrum  must  be 
obtained  over  one  short  frequency  interval  at  a  time,  and  cannot 
usually  be  displayed  in  a  single  operation  as  is  the  case  in 
high  resolution  nuclear  magnetic  resonance.  The  other  two 
advantages,  however,  have  provided  the  impetus  for  our  development 
of  FT-NQR  techniques. 

Once  FID  or  echo  signals  have  been  obtained,  the  proper  FT 
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treatment  for  each  will  yield  the  desired  lineshape  spectrum  ' 

In  this  section  the  procedures  the  authors  found  most  satisfactory 
are  discussed. 
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The  problem  of  phase  correction  for  NQR  spectra  presents  a 
particular  problem  not  found  in  f ix?d-frequency  spectrometers.  By 
phase  correction  here  we  do  not  mean  the  usual  procedure  performed 
on  FID  data  to  correct  for  the  first  few  missing  data  points  lost 
in  the  spectrometer  recovery  time.  Rather,  we  are  concerned  with 
the  problem  present  in  both  FID  and  echo  data  that  the  output  of 
our  phase-sensitive  receiver  is  a  generally  unknown  admixture  of 
"absorption"  and  "dispersion"  signal.  The  degree  of  admixture 
varies  in  a  complicated  manner  with  the  spectrometer  operating 
frequency  and  tuning.  The  phase  shift  of  the  phase-sensitive 
receiver  output  away  from  true  "absorption"  can  be  considerable 
and  must  be  mathematically  corrected  in  order  to  obtain  a  true 
line  shape  from  the  time-domain  data.  In  Table  2.2-1  the 
mathematical  results  for  Lorentzian  lineshapes  are  collected  as 
guides  to  the  solution  of  this  problem.  (A  more  elegant  solution 
would  be  possible  if  quadrature  detection  were  available.) 

The  results  in  the  table  show  that  the  presence  of  phase 
shifts  produces  an  admixture  of  absorption  and  dispersion  modes  in 
the  case  of  FID  signals,  and  a  possible  loss  of  intensity  in  echo 
signals.  This  problem  can  be  avoided  while  at  the  same  time  the 
true  Lorentzian  line  shape  is  preserved,  if  the  modulus  squared 
transform  is  computed  for  FID  signals,  and  the  modulus  transform 
is  computed  for  echoes.  Thus  the  proper  line  shapes  are  obtained 
in  each  case  regardless  of  the  degree  and  source  of  phase  shift 
without  need  for  a  separate  "phase  correction"  subroutine. 

In  Figures  2.2-1  to  2.2-8  the  results  of  computer-simulated 
soectra  are  known.  An  ecno  and  an  FID  signal  have  been  simulated 
for  both  Lorentzian  and  Gaussian  lineshapes.  Cosine,  Sine,  the 
square  of  the  modulus,  and  modulus  transforms  are  computed  and 
displayed.  Figures  2.2  1  to  2.2-4  are  for  zero  phase  shift,  while 
Figures  2.2-5  to  2.2-8  have  a  phase  shift  of  30°  in  the  time 
domain  signals.  Note  that  in  all  cases  the  conclusions  discussed 
for  proper  data  processing  are  borne  out,  mainly  that  the  modulus 
squared  transform  should  be  used  for  FID  signals,  and  the  modulus 
transform  should  be  used  for  echo  signals. 
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TABLE  2.2-1 


LORENTZIAN  FID  AND  ECHO  COMPLEX 
SIGNALS  AND  THEIR  FOURIER  TRANSFORM 


FID 


TIME  DOMAIN:  f(t)  =  e^W  e"1<t>  t  >  0 


C,  Cosine  transforra  =  — - - 2  |acos<j>  -  (aj-coQ ) si n4> } 


a  +(w-q>0) 


S,  Sine  transform  =  —g — - 2  { ots i n ^  +  (uj-oo0 )cos<f>  | 


2  2 

C  +S  ,  Modulus  squared  transform 


2  2  * 

[C  +S  ]  ,  Modulus  transform 


ct  +(u-u iQ)‘ 
1 

a2+(u)-m0)2 
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ECHO 
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e'ate  °Jot  e'1* 


C,  Cosine  transform  * 


a2+{w-aj„)2 
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S,  Sine  transform 
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2  2 

C  +S  ,  Modulus  squared  transform 


[C2+S2]*,  Modulus  transform 
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Figure  2.2-1 


Computer-simulated  exponential  FID  signal 
(a)  and  its  Fourier  transform  (b-e) .  (a) 

time  domain  signal  [e-atcos  wt] ;  (b)  Cosine 
FT0°f  (3j;  (c)  Sine  FT,  (d)  modulus  squared 
(C2  +  S<=0  (e)  modulus  FT,  /c2  +  s2‘ 


Fi«4*4*CV  (Hi! 

Id )  Uodutui  Sfluortd  FT 


Ft*e«««ct  l"ti 
(elMoCului  FT 


re  2.2-2.  Computer-simulated  exponential  echo 

signal  (a)  and  its  Fourier  transforms 
(b-e| .  (a)  time  domain  signal 

[e-alt  I  cos  ut) ;  (b)  Cosine  FT  of  (a) ; 
(c)  Sine  FT  vanishes  in  this  case;  (d) 
modulus  squared  FT,  (c2  +  s2) . 

(e)  modulus  FT,  /c2 
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(d)  Uodabil  SQMftd  FT 


(Hji 

(*)  UoduUit  FT 


2.2-3.  Computer-simulated  Gaussian  FID  signal 
(a)  and  its  Fourier  transforms^ (b-e) . 

(a)  time  domain  signal  “cos  ut^ti-Ol 

(b)  Cosine  FT  of  (a) ;  (c)  Sine  FT, 

(d)  modulus  squared FT, (C2  +  s2) 

(e)  modulus  FT,  /c~2  +  s2 
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Figure  2.2 
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4.  Computer- simulated  Gaussian  echo  signal 
(a)  and  its  Fourier  transforms  (b-e) . 

(a)  time  domain  signal  (e~(at)2  cos  wt]  ; 

(b)  Cosine  FT?  (cj  Sine  FT  vanishes  in 
this  case;  (d)  modulus  squared  FT,  (C2+  s?) 
(e)  modulus  FT  ,  »'c2  +  s2 
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Figure  2.2-5.  Computer- simulated  exponential  FID  signal 
with  phase  shift  (a)  and  its  Fourier 
transforms  (b-e) .  (a)  time  domain  signal 

(e-atcos  <wt  -  $);  (b)  Cosine  FT  of  (a) ; 

(c)  Sine  FT;  (d)  modulus  squared  ft,  (C2+  £2) 
(e)  modulus  FT,  Jc2  '+  s2 
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ure  2.2-6.  Computer-simulated  exponential  echo  signals 
with  phase  shift  (a)  and  its  fourier  trans¬ 
forms  (b-e) .  (a)  time  domain  signal 

fe-altfcos  (wt  -  6)1;  (b)  Cosine  FT  of 
(a);  \c)  Sine  FT;  (d? modulus  squared  FT, 
{C2*  $2)  (e)  modulus  FT,  »'&+  S* 


Figures  2.2-9,  2.2-10  and  2.2-11  show  experimental  spectra 
that  demonstrate  the  foregoing  arguments.  Figures  2.2-9 (a)  and 

2.2- 10 (a)  are  the  Nitrogen-14  NQR  FID  signals  at  77°K  from 
hexamethylenetetramine  (HMT)  and  urea,  respectively.  The  Cosine 
and  Sine  transforms  of  HMT,  Figures  2.2-9 (b)  and  (c) ,  clearly  show 
the  admixture  of  absorption  and  dispersion  expected  when  phase 
shifts  exist  in  the  time-domain  data.  This  effect  is  much  less 
evident  in  the  cosine  and  sine  transforms  of  urea,  Figures 

2.2- 10 (b)  and  (c) ,  which  occured  with  only  a  small  phase  shift. 
Finally,  the  proper  lineshapes  are  shown  in  Figures  2.2-9 (d)  and 

2 . 2- 10 (d) ,  the  modulus  squared  transforms  of  the  time-domain 
signals.  Note  the  fine  structure  on  the  HMT  line,  first  reported 
by  Colligiani  and  Ambrosetti.  The  modulus  transforms,  Figures 

2.2- 9  (e)  and  2.2-10  (e) ,  are  shown  for  comparison  and  they  are 
visibly  broader  than  the  true  lineshapes. 

An  example  of  an  echo  signal  is  shown  in  Figure  2.2-11,  a 
doublet  of  v_  lines  from  the  monoclinic  phase  of  TNT  at  77°K. 
Figures  2.2-JKb)  and  (c)  are,  respectively,  the  Cosine  and  Sine 
transforms  of  this  signal.  Note  that,  as  expected,  both  of  these 
transforms  yield  valid  lineshapes  with  fractional  amplitudes,  and 
that  the  modulus  transform,  Figure  2.2-11  (e)  ,  yields  the  correct 
lineshape  with  maximum  amplitude. 

In  Appendix  B  we  collect  a  library  of  TNT  lineshapes.  In 
each  case  the  time-domain  echo  signal  is  shown  at  the  top  of  each 
figure,  and  the  modulus  Fourier  transform,  i.e.»  the  lineshape,  is 
shown  at  the  bottom  of  the  figure. 
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Figure  2.2-9  Experimental  N  FID  signal  from 
hexamethylenetetramine  (HMT)  at 
room  temperature  (a)  and  its  Fourier 
transforms  (b  -  e)  .  (a)  time  domain 

sictnal,  (b)  Cosine  FT  of  fa);  (c) 

Sine  FT;  (d)  modulus  squared  (c2+s2) 
(e)  modulus  FT,  /c'2+s* 
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Figure  2.2-10,  Experimental  H14  FID  signal  from 
the  v+  line  in  urea  at  77’K  a  and 
its  Fourier  transforms  (b-e) .  (a) 

time  domain  signal,  b)  Cosine  FT  of 
(a) j  (c)  Sine  FT?  (d)  modulus  squared  FT, 
(c2+s2)  (e)  modulus  FT,  /c2+s2 


2-31 


Miihttconds 
(o)  Time  Domain 


Figure  2.2-11  Experimental  N  spin  echo  signal  from 

the  727  kHz  v_  doublet  in  monoclinic  TNT 
at  77°K  (a)  and  its  Fourier  transforms 
(b  -  e) .  (a)  time  domain  signal,  (b)  Cosine 
FT  of  (a);  (c)  Sine  FT;  (d)  modulus  squared  FT, 
(C2+S2)  {e)  modulus  FT,  /c2+S2 


In  the  last  section  it  was  pointed  out  that  unlike  the  case 
of  high-resolution  nuclear  magnetic  resonance  (NMR)  it  is  not 
generally  possible  for  a  single  Fourier  transform  (FT)  to  display 
a  complete  NQR  Fourier  transform  (FT) .  In  fact,  typical 
Nitrogen-14  NQR  Fourier  transforms  (FT)  span  only  a  few  kHz.  For 
examples,  see  Appendix  B.  There  are  several  reasons  for  this 
restriction: 

A.  RECEIVER  BANDWIDTH 

In  order  to  maximize  the  signal-to-noise  ratio  (S/N) ,  the 
receiver  bandwidth  is  kept  low,  typically  not  much  larger  than  10 
kHz.  In  any  event,  the  cut-off  frequency  must  not  be  higher  than 
half  the  frequency  of  the  analog-to-digital  converter  of  the 
signal  averager  used. 

B.  SAMPLE  COIL  QUALITY  FACTOR 

High  values  of  the  sample  coil  quality  factor  ,  Q,  are  used  to 
maximize  the  S/N.  The  half-width  at  half- maximum  of  the 
receiver  coil  is  given  by  Af  =  fQ/Q.  Thus,  typical  values  of  the 
transmitter  frequency  f  =  1MHz  and  Q  =  100  yield  Af  =  10  kHz. 

Thus  we  see  that  under  these  conditions  the  sample  coil,  quality 
factor  introduces  considerable  distortions  for  Fourier  transforms 
that  extend  much  beyond  20  kHz. 
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C.  TRANSMITTER  PULSE  SHAPE  FUNCTIONS 


The  spectral  intensity  of  rectangular  transmitter  pulses  of 
duration ,  t,  and  carrier  RF  frequency,  fQ  =  uQ/2-n ,  varies  with 
frequency  as  (rAft)  sin  (rAft)  ,  where  Af  is  the  frequency  offset 
from  fQ.  Thus,  there  is  no  transmitter  power  whatever  at  the 
first  zero  of  this  function,  i.e.,  at  a  frequency  of  Af  =  1/t  away 
from  the  carrier  fQ.  Thus,  the  width  of  an  NQR  spectrum  ,D ,  that 
is  to  be  evenly  excited  by  one  or  more  transmitter  pulses  should 
be  small  compared  to  1/t.  Remembering  that  the  choice  of  pulse 
width  is  governed  by  the  need  to  generate  an  optimum  flip-angle 
(cf.  Equation  2.1-1),  these  considerations  place  a  rigid 
constraint  on  D.  For  our  spectrometer5,  which  is  comparable  to 
other  designs  in  use  in  this  respect,  t  s  50  microseconds.  Thus, 
one  requires 


D  <<  20  kHz, 

which  is  clearly  a  very  stringent  requirement. 


(2. 2. 1-1) 


The  above  considerations,  together  with  the  fact  that  the 
typical  NQR  spectrum  reported  in  the  literature  consists  of  only  a 
few  widely  separated  narrow  lines,  have  resulted  in  FT/ NQR 
spectral  widths, D,  of  the  order  of  3  to  5  kHz.  However,  driven  by 
a  desire  to  display  as  much  as  possible  of  the  crowded  14N  NQR 
spectrum  of  TNT  (cf.  Figure  3.1-1  in  the  next  section,)  we  were 
led  to  a  study  of  ways  to  broaden  D.  The  following  steps  were 
taken, 


a.  The  signal  average  analog-to-digital  (A-D)  converter 
speed  was  increased  from  50  kHz  to  200  kHz.  This  was 
accomplished  by  using  a  Princeton  Applied  Research  PAR 
4202  signal  averager  with  a  9  bit,  200  kHz  A-D  converter 
and  by  designing  and  building  a  custom  interface  for  the 
Digilab  Data  System  NOVA  minicomputer.  (Details  of  this 
work  are  given  in  Appendix  d ) . 
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b.  The  spectrometer  receiver  output  filters  were  changed 
to  allow  the  choice  of  a  100  kHz  low  pass  filter. 

c.  The  sample  coil  Q  was  lowered  to  a  value  of  24  by  adding 
a  series  resistance  to  it. 

d.  Transmitter  pulse  widths  were  reduced  to  10 
microseconds.  Because  of  finite  rise  and  fall  times, 
this  value  could  not  be  reduced  further. 

These  modifications  allow  spectra  with  width  D  z  50  KHz  in  a  single 
FT  to  be  displayed,  roughly  an  order  of  magnitude  improvement  over 
past  practice.  Signal  averaging  times  needed  to  be  sharply 
increased  to  make  up  for  the  degraded  signal  to  noise  ratios 
resulting  from  the  choices  made  in  "a"-  "d"  above.  In  order  to 
minimize  distortion  of  the  computed  broadband  spectrum  the 
transmitter  pulse  RF  carrier,  f  ,  was  placed  at  the  edge  of  the  ] 

desired  spectrum,  D,  while  the  receiver  and  tne  sample  tank  circuit  ; 

were  tuned  to  the  middle  of  D.  The  spectrometer  was  designed5  to  j 

have  a  low  value  of  Q,  about  10,  during  the  transmitter  pulse.  < 

Therefore,  placing  fQ  at  the  edge  of  D,  the  position  of  choice  to  J 

obtain  the  true  spectrum  in  the  absence  of  quadrature  detection 
poses  no  significant  danger  of  distortion  due  to  the  sample  coil 
Q. 

Figure  2. 2. 1-1  shows  our  result  for  monoclinic  TNT  at  77°K. 

All  six  v+  lines  are  contained  in  a  single  interferogram,  which  j 

is  itself  shown  at  the  top  of  the  figure.  The  final  signal  was  j 

obtained  by  the  coherent  addition  of  600,000  echoes  as  follows:  a 
SLSE  train  was  generated  with  100  echoes  at  10  millisecond 
separation  each.  The  spin  system  was  allowed  to  approach  thermal 
equilibrium  with  the  lattice  for  3  seconds,  and  the  process 
repeated  6,000  times.  The  total  observation  time  was  6.7 
hours.  The  transmitter  frequency  and,  therefore,  also  the  origin 
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Figure  2.2. 1-1  Echo  signal  interferogram  and  its  broadband 

Fourier  transform  showing  all  oix  v+  lines  of 
monoclinic  TNT  at  77°K.  Other  parameters  are 
given  in  the  text. 


Various  corrections  can  be  made  to  the  raw  data  shown  in 

Figure  2. 2. 1-1  in  an  attempt  to  improve  the  considerable  distortion 
still  seen  there,  exemplified  by  the  decreasing  amplitude  of  peaks 
with  increasing  frequency  offset.  We  list  some  possible  remedies, 


a.  correct  for  sample  coil  Q  function 

b.  correct  for  the  receiver  filter  characteristics 

c.  correct  for  the  transmitter  pulse  spectral  intensity 
variations 

d.  correct  for  the  transmitter  pulse  non-ideality,  e.g., 
the  phase  shifts  introduced  by  finite  rise  and  fall 
times. 

e.  correct  for  the  fact  that  the  intensity  of  an 
interferogram  signal  obtained  from  the  coherent  sum  of  a 
train  of  echoes  ir.  a  SLSE  sequence  is  a  function  of  the 
relaxation  times  of  the  individual  resonances  and  not 
merely  a  reflection  of  a  density-of-states  function. 


In  Figure  2.2. 1-2  we  show  the  same  TNT  spectrum  of  the  previous 
figure  corrected  only  for  a.  and  b.  above.  A  slight  improvement 
is  observed.  (N.B.,  the  sharp  feature  exactly  at  50.0  kHz  is  a 
spurious  "electronic  resonance"  and  is  not  sample  dependent) .  The 
correction  function  was  obtained  experimentally  by  measuring  the 
response  of  the  coil  and  receiver  combination  in  a  separate 
procedure.  The  resulting  normalizing  function  is  shown  as  Figure 
2.2. 1-3. 
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Figure  2. 2, 1-3  Normalizing  function  for  Broadband 

FT  of  TNT  which  was  used  to  obtain 
Figure  2.2-2  from  Figure  2.2-1. 
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One  of  the  benefits  of  being  able  to  display  a  larger 
spectrum  D,  is  that  it  enables  us  to  study  effects  that  occur  in 
this  larger  domain.  An  obvious  case  is  the  situation  of  the  need 
to  study  NQR  lines  that  are  themselves  intrinsically  broad  in 
frequency  space.  An  example  is  the  hydrogen  bonded  crystal 
Hexamethylenetetramine  Hexahydrate  which  has  three  NQR  lines 
labelled  v+,  v_,  and  vQ  broaded  by  positionally  disordered  protons 
in  the  structure.  The  broadband  FT  spectra  of  these  lines  are 
shown  in  Appendix  C . 

2 .3  TEMPERATURE  DEPENDENCE 

Two  methods  were  used  to  determine  the  frequency/temperature 
dependence  of  the  TNT  NQR  lines. 

The  first  method  of  temperature  control  is  shown  in  Figure 
2.3 -la,  the  temperature  stabilization  assembly.  This  assembly 
with  the  circuitry  shown  in  Figure  2.3-2,  was  able  to  control  the 

temperature  to±0.1°C  between  140°K  and  320°K.  For  temperatures 
reasonably  above  ambient »  the  liquid  nitrogen  ballast  is  not 
needed . 
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DESCRIPTION _  j  DATE  1  APPROVED 


The  temperature  control  circuitry  proportionally  controls, 
and  turns  "on"  only  during  a  zero  crossing  of  the  AC  signal  to 
minimize  RF  interference.*-  No  crosstalk  into  the  instrument  was 
observed. 

The  Nalgene® pipet  wasner  did  fairly  well  holding  liquid 
nitrogen.  It  is  necessary  to  fill  it  somewhat  slowly  in  order  to 
minimize  the  thermal  shock. 

A  second,  cruder,  method  was  also  used  in  the  initial  runs. 
This  was  a  brute  force  approach  of  cooling  the  sample  to  77 °K  via 
liquid  nitrogen,  and  following  the  frequencies  as  the  sample 
slowly  warmed  to  room  temperature  in  a  Dewar.  The  apparatus  used 
is  shown  in  Figure  2.3-lb,  the  temperature  run  assembly.  Also 
shown  is  a  typical  Temperature  vs.  time  plot. 


Applications  of  the  yA742  TRIGAC ,  A  Zero  Crossing  AC  Trigger 
by  R.B.  Hood,  1970  Fairchild  Semiconductor  04-10-0063-50. 

2nd  printing. 


©Registered  trademark  of  the  NALGE  Corp. 
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2.4  A  METHOD  TO  PAIR  UP  V_  AND  v+  LINES 

In  this  section  we  outline  a  simple  double  resonance  method 

of  pairing  up  v+  and  v_  lines  corresponding  to  the  same  NQR  site. 

A  uniform  static  magnetic  field  Hq  is  applied  to  the  sample  and 

a  spin-locked  spin  echo  (SLSE)  experiment  is  performed.  The 

spacing  between  pulses  is  chosen  small  enough  that  the  echo 

train  decay  constant  T2c  is  appreciably  longer  than  T2*  The 

individual  echoes  in  the  echo  train  are  co-added  to  give  a 

resultant  signal  whose  magnitude  is  then  plotted  vs.  the 

magnetic  field  Ho.  For  a  unique  value  of  Hq,  the  NMR  splitting 

1 4 

of  protons  m  the  neighborhood  of  the  N  site  exactly  matches 
14 

the  N  "difference"  line  vq=v+-v_.  When  this  happens,  mutual 

spin  flips  between  the  Proton  and  Nitrogen  spin  systems  will 

14 

result  in  lower  values  of  the  N  relaxation  times  T2  and  T^ . 
Since  the  inequality  T1>T2e>T2,  holds,  the  concomitant 
decrease  m  Tn  results  in  an  often  dramatic  decrease  in 
the  value  of  the  resultant  co-added  echo. 

Figure  2.4-1  shows  just  such  a  plot  of  NQR  echo  height 
vs.  the  proton  NMR  frequency  yHQ,  where  y  is  the  proton 
magnetogyric  ratio.  The  NQR  transition  under  study  is  the 
811  kHz  v_  lines  in  m-dinitrobenzene.  The  sharp  dip  in  the 
echo  enabled  the  paring-up  of  this  line  with  the  previously 
detected  v+  line  at  1030  kHz. 
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Figure  2.4-1. 


FREQUENCY 

Plot  of  the  amplitude  resultant  from  coherently 
adding  echoes  in  a  single  SLSE  train  vs.  the 
applied  static  magnetic  field  expressed  as  the 
corresponding  proton  NMR  resonance  frequency. 
The  NQR  transition  is  the  811  kHz  v_  line  of 
m-dinitrobenzene  at  77°K. 
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3 . 0  NQR  STUDY  OF  TNT 


It  has  long  beer,  established  that  a-Trinitrotoluene  can 
19  2  0 

coexist  '  as  two  different  phases  at  ambient  temperatures.  A 

monoclinic  form  is  the  product  of  annealing  cast  TNT  and  has  been 

212  3 

studied  by  several  workers,  '  A  crystal  structure  has  only 

24 

recently  been  obtained  by  J.R.C.  Duke.  Another  possible  crystal 

232527 

habit  has  orthorhombic  symmetry  '  '  .  This  phase  can  be  formed 

by  quickly  chilling  molten  TNT  as  well  as  by  recrystallization 

from  cyclohexanol.  We  have  recently  become  aware  that  a  crystal 

28 

structure  of  this  phase  also  has  been  obtained  by  J.R.C.  Duke 

Samples  of  cast  TNT,  which  commonly  contain  a  mixture  of 
these  two  phases,  will  often  develop  cracks  and  voids  when 
orthorhombic  microcrystals  transform  to  the  denser  monoclinic 
phase  upon  annealing.  The  desire  to  learn  more  about  the 
molecular  structure  and  molecular  dynamics  of  TNT  in  the  two 
phases  provided  part  of  the  rationale  for  this  study. 
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3.1  SAMPLE  PREPARATION  AND  NQR  SPECTRUM 


a-TNT  ootained  from  Eastman  Kodak  arrives  in  the  form  of 
small  flakes  moist  with  10%  water.  These  flakes  were  gently  dried 
under  a  heat  lamp  and  recrystallized  from  acetone.  The  resulting 
very  pale  yellow  crystals  were  used  to  prepare  samples  of  solid 
TNT  which  were  preferentially  crystallized  in  one  of  its  two 
crystal  habits.  Monoclinic  TNT  was  obtained  by  seeding  the 
melt^'^O/^l  with  about  h  of  a  percent  of  2 , 2 1  , 4 ,4  '  , 6 , 6 ' - 
-Hexanitrostylbene  (HNS)  allowing  the  complex  (TNT)2*HNS  to  form 
and  then  to  nucleate  the  growth  of  monoclinic  TNT. 

The  resulting  polycrystalline  material  consisted  of  a  mass  of 
unoriented  and  very  small  crystallites  as  expected.  Orthorhombic 
TNT  was  prepared  by  quenching  molten  TNT  in  small  batches  on  a 
copper  slab  maintained  at  about  0°C. 
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The  N  NQR  spectrum  of  monoclxnic  TNT  and  of  orthorhombic 

TNT  is  shown  in  Figure  3.1-1.  The  spectrum  of  each  phase  consists 

of  twelve  resonance  lines,  corresponding  to  six  distinct  nitrogen 

sites.  Each  site  gives  rise  to  a  v_  and  a  corresponding  v+  line 

(See  Appendix  A)  which  are  related  to  the  quadrupole  coupling 
2 

constant,  e  qQ/h  ,  and  to  the  asymmetry  parameter, n,  by  the 
following  equations: 

e^qQ/h  =  -  (v  +  v  ) 

(3 . 0-] 

(v  _  v  ) 
n  -  3  — + - 3- 

(V+  +  V_) 
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For  each  corresponding  pair  of  v_  and  v+  lines  there  also 
exists  a  "difference  line"  v ^  whose  value  is  given  by  =  v+-  v_. 

In  TNT  the  frequencies  of  these  transitions  fall  below  the  lower 

operating  point  of  our  spectrometer  and  were  not  detected. 

2 

Therefore,  in  order  to  compute  e  qQ/h  and  h  for  each  site  it  is 
necessary  to  "match-up"  the  corresponding  v+  and  v_  lines.  This 
was  done  using  a  simple  double  resonance  technique  described  in 
Section  2.4.  Table  3.1-1  shows  the  NQR  frequencies  matched  and 
assigned  to  six  crystallographic  sites  for  both  phases  of  TNT. 

The  similarity  between  the  NQR  spectra  of  the  two  phases  of 
TNT  as  depicted  in  Figure  3.1-1  is  remarkable.  Experience  with  NQR 
spectra  of  other  subtances  which  undergo  structural  phase  transition 

shows  that  the  NQR  frequencies  undergo  shifts  of  the  order  of  10%. 
Now,  these  structural  phase  transition  shifts  have  been 
interpreted  as  giving  a  measure  of  the  contribution  of  atoms  and 
molecules  other  than  those  directly  bonded  to  the  resonant  site  to 
the  total  field  gradient  at  that  site.  An  independent  measurement 
of  the  same  effect  is  the  size  of  "crystal  splitting",  the 
difference  in  NQR  frequencies  corresponding  to  sites  that  are 
chemically  equivalent,  but  which  are  made  physically  inequivalent 
by  the  crystal  packing.  These  two  measures  of  crystal  field 
effects  have  been  found  generally  in  good  agreement  when  they 
could  be  compared. 

In  the  present  case,  we  observe  that  the  spectra  of 
orthorhombic  and  monoclinic  TNT  are  nearly  identical  to  one 
another  while  at  the  same  time  each  phase  displays  considerable 
physical  inequivalence  (See  Section  3-2) .  This  allows  us  to 
conclude  with  some  confidence  that  the  structure  and  conformation 
of  each  molecule  in  orthorhombic  TNT  and  its  relation  to  its 
immediate  molecular  neighbors  must  be  remarkably  similar  to 
the  case  in  monoclinic  TNT. 
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TABLE  3.1-1 


Nitrogen-14  NQR  Spectra  of  Trinitrotoluene 
at  77°K.  The  spectral  lines  are  assigned 
to  one  para  (p)  site  and  two  ortho  sites 
(oA,oB)  for  each  of  two  inequivalent  mole¬ 
cules  (I  and  II).  Data  for  the  monoclinic 
phase  are  listed  first,  with  corresponding 
lines  of  the  orthorhombic  phase  shown  in 
parentheses. 


MOLECULAR  SITE 

V 

kHz 

v  + 
kHz 

oAI 

792.1 

(792.4) 

869.4 

(867.8) 

oAII 

801.8 

(803.2) 

895.4 

(897.9) 

oBI 

767.3 

(766.5) 

875.6 

(875.7) 

oBII 

767.7 

(768.5) 

857.0 

(858.6) 

Pi 

729.7 

(727.4) 

861.8 

(862.3) 

pll 

_ 

730.1 

(732.5) 

888.1 

(885.6) 
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Because  of  this  great  similarity,  our  discussion  of  structure 
and  bonding  of  ct-TNT  will  focus  specifically  on  the  temperature 
stable  (monoclinic)  phase  of  this  material. 

3.2  Structure  and  Bonding  of  TOT  -  General 

For  each  crystal  phase  of  TNT  the  observation  of  six  v  lines 
and  six  v+  lines,  indicating  the  presence  of  six  Nitrogen  NQR 
sites  is  in  agreement  with  the  x-ray  crystallography  work  of  J.R.C. 
Duke2^'2^.  There  are  two  crystallographically  distinct  molecules 
per  unit  cell,  each  with  physically  inequivalent  nitro-groups . 
Figures  3.2-1  and  3.2-2  shows  Duke's  results  for  the  molecular  and 
crystal  structure  for  the  two  phases  of  TNT. 

The  isolated  molecule  of  trinitrotoluene  (2-methyl-l , 3,5- 
trinitrobenezene)  has  two  chemically  distinct  nitro-groups  which 
we  could  designate  p  and  o  for  the  groups  which  are  para  and 
ortho,  respectively,  to  the  methyl  group.  If  the  two  ortho  sites 
become  physically  inequivalent  in  the  solid,  and  we  label  them  oA 
and  oB  respectively,  and  if  we  now  label  the  two  crystallo¬ 
graphically  distinct  molecules  per  unit  cell  I  and  II,  we  reach  a 
nomenclature  for  the  six  sites  which  we  will  use  for  the  remainder 
of  this  report.  Namely,  pi,  oAI  and  oBI  for  molecule  I,  and  pH, 
oAII  and  oBII  for  molecule  II. 
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Figure  3.2-1.  Projection  of  the  crystal  and  molecular 
structure  of  MONOCLINIC  TNT  on  the  a-c 
crystal  plane.  Data  and  figure  obtained 

from  J.R.C.  Duke  (private  communication, 
1982)  . 
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Figure  3.2-2.  Projection  of  the  crystal  and  molecular 
structure  of  ORTHORHOMBIC  TNT  on  b-c 
crystal  plane.  Data  and  figure  obtained 
by  J.R.C.  Duke  (private  communication, 
1982) 


BLOCK 


3-8 


The  above  considerations  about  the  chemical  and  physical 
(in) equivalence  of  NQR  sites  is  usually  extremely  useful  in 
comparing  these  sites  to  their  corresponding  crystallographic 
sites.  However,  this  is  not  the  case  in  TNT  as  can  be  seen  by  the 
seemingly  incongruous  pairing  of  and  v+  lines  in  Figure  3.1-1. 

In  fact,  it  will  be  shown  that  the  chemical  inequivalence  of  "0"  and 
"p"  sites  are  not  basically  relevant  to  tne  interpretation  of  NQR 
spectrum  of  Figure  3.1-1.  Rather,  the  determining  facte  in 
explaining  the  splitting  of  the  NQR  lines  is  the  physical 
inequivalence  of  the  six  nitrogen  sites  in  the  solid  state. 
Specifically,  we  will  show  that  the  parameter  that  makes  the 
dominant  contribution  is  the  angle  of  twist  of  the  nitro-group 
about  the  C-N  bond. 

In  the  following  sections,  we  will  first  develop  a  set  of 
equations  (Townes-Dailey  method)  relating  the  observed  NQR 
parameters  to  a  bonding  model  of  the  nitro-group.  Then  we  will 
use  these  equations  to  interpret  the  TNT  data,  and  show  how 
additional  NQR  data,  much  of  it  new,  on  molecules  similar  to  TNT 
supports  our  interpretation. 

3.2.1  Structure  and  Bonding  of  TNT: Townes-Dailey  Theory 

32 

The  Townes-Dailey  (TD)  method  of  relating  measured  field 
gradients  to  orbital  occupation  numbers  is  a  first-order  theory 
which  works  best  when  comparisons  are  made  among  similar 
molecules33.  In  our  case,  where  we  limit  our  discussion  to 
C-nitro  groups,  we  expect  that  the  analysis  will  provide  useful 
insights  in  the  variation  of  structure  and  electron  distribution 
among  the  various  nitro  groups  under  study. 
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Figure  3. 2. 1-1  depicts  the  coordinate  system  which  is 

used  to  analyze  the  bonding  of  the  nitro  group.  Because  two  of 

the  axes  are  perpendicular  to  local  mirror  planes,  one  is  assured 

that  the  chosen  axes  will  diagonalize  the  electric  field  gradient 

due  to  the  local  electron  distribution.  The  axes  are  labelled 

using  the  convention  independently  arrived  at  by  Subbarao  and 
34  35 

Bray  ,  and  by  Cheng  and  Brown  .  The  nitrogen  site  is  assumed 

**\ 

to  be  in  a  state  of  sp*"  hybridization:  a  sigma  bonding  orbital  of 

occupation  number  a forms  a  bond  to  the  ring  carbon  atom,  while 

the  two  remaining  hybrids,  each  with  occupation  number  owQ,  form 

sigma  bonds  to  the  oxygen  atoms.  A  nitrogen  orbital,  with 

occupation  number  ir  ,  takes  part  in  the  molecular  pn  orbitals. 

Table  3. 2.1-1  lists  the  nitrogen  hybrid  orbitals  and  their 

2 

occupation  numbers.  The  standard  sp  hybridization  is  used  in 
these  equations  even  though  the  angle  0N0  is  known  to  be  larger 
than  120°  by  about  four  to  five  degrees.  Although  the  ,?D  method 
can  easily  accommodate  this  "refinement"  it  has  been  the 
experience  of  this  and  other  workers  that  it  does  not  result  in  a 
significant  improvement  of  the  final  results. 


Using  the  orbitals  of  Table  3. 2.1-1,  the  Townes-Dailey 
relations  can  be  written  in  the  following  very  useful  form. 


V%  =  3  <*  *  “no’ 


(3.2-la) 


v-/qo  =  2  *°NC  "  °NO) 
or  equivalently  as. 


(3.2-lb) 


t’-'W  =  *  v+/q' 


{oNC  ~3N0)  =  2  v~/qo 


(tJNC  ■*)  =  (2v“  '  t  v+)/q 


(3.2-2a) 

(3 . 2-2b) 

(3.2-20 
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TABLE  3.2. 1-1 


2 

sp  hybrid  orbitals  and  occupation  numbers 
for  the  nitrogen  site  in  C-N i tro compounds. 
The  coordinate  system  of'  Figure  3. 2. 1-1 
is  used. 


— 

OCCUPATION  NUMBER 

NITROGEN 

ORBITAL 

°N0 

^NOI 

■/Js* 

4"*  + 

4p2 

aN0 

%)II 

=.4s  + 

4px- 

4pz 

.  aNC 

*NC 

=  4s  ~ 

/Ip 

v  3 

TT 

=  py 

~ — — _ 
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Here  qQ  is  the  magnitude  of  the  coupling  constant  per  p-electron, 
a  parameter  estimated  to  be  between  -8  and  -10  MHz.  Although  most 
of  our  conclusions  will  be  insensitive  to  the  particular  choice 
of  q  ,  we  prefer  the  latter  value,  and  shall  use  it  when  necessary. 
The  sign  of  qQ  has  been  chosen  so  that  all  forms  in  Equations 
3.2-2  are  positive. 

3.2.2  Structure  and  Bonding  of  TNT:  Nitro  Group  Conformation 

Equations  3.2-2  are  the  basis  for  the  interpretation  of 
the  NQR  data  of  ct-TNT  in  terms  of  the  structure  and  bonding  of 
this  molecule. 

The  first  effect  to  be  considered  in  the  differential  effect 

which  the  methyl-group  will  have  on  nitro-groups  ortho  and  para  to 

it.  Fortunately,  this  general  question  of  substituent  effects  on 

NQR  parameters  has  been  thoroughly  investigated,  most  notably  by 

33 

Bray  and  coworkers  and  is  well  understood.  The  main  results  is 
that  good  correlations  can  be  obtained  between  measured  field 
gradients  and  Hamroet  a  constants.  The  later  are  a  measure  of  the 
ability  of  the  substituent  to  withdraw  or  supply  electrons  to  a 
reaction  site.  It  should  be  noted,  however,  that  these  constants 
are  determined  in  the  liquid  state.  Thus  it  is  not  surprising  that 
their  correlation  with  NQR  data  obtained  in  solids  breaks  down 
when  anisotropic  interactions  are  strong. 

In  the  case  of  a~TNT,  considerations  of  the  chemical 
inequivalence  of  the  ortho  and  para  nitro-groups  using  the 
correlcations  found  by  Subbarao  and  Bray  ,  lead  to  expectations  of 
fx*equency  shift  less  than  10  kHz.  Clearly,  these  shifts  are 
much  less  than  the  splittings  of  actual  TNT  frequencies.  We 
conclude  that  chemical  inequivalence  due  to  the  electron  releasing 
property  of  the  methyl  group  cannot  even  begin  to  explain  the 
observed  NQR  data  in  &-TNT. 


SLOCK 
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We  now  turn  to  a  consideration  of  the  possibility  of  strong 
physical  inequivalence  of  the  nitro-groups  in  nitrobenzenes .  In 
fact  crystallographic  studies  of  these  substances  show  that  the 
plane  of  the  N02~group  can  exist  in  equilibrium  in  solids  at 
relatively  large  angles  of  rotation  from  the  plane  of  the  benzene 
ring.  This  rotation  occurs  mainly  about  the  CN  bond,  not 
withstanding  the  large  expected  double  bond  character  of  this 
bond.  The  crystallographic  studies,  taken  together,  indicate  that 
this  angle  of  twist  generally  increases  with  the  number  of 
neighboring  substituents  on  the  benzene  ring.  What  is  surprising 
is  that  even  with  no  neighbors,  angles  of  up  to  10°  or  more  are 
common.  One  substituent  ortho  to  the  nitro-group  can  lead  to 
twist  angles  of  the  order  of  30°  to  40° ,  while  rotations  of  up  to 
60°  or  70°  are  found  when  both  sites  ortho  to  the  nitro-group  have 
a  bulky  substituent.  In  addition  to  the  effect  of  steric 
hindrance  from  neighboring  substituents,  the  details  of  crystal 
packing  and  of  intermolecular  interactions  are  responsible  for 
further  variations  in  the  nitro-group  twist  angles. 

When  the  plane  of  a  nitro-group  in  a  nitro  benzene  is 
constrained  to  be  twisted  away  from  the  plane  of  the  benezene  ring 
by  a  rotation  about  the  CN  (for  any  of  the  reasons  outlined  above) 
it  follows  that  at  the  same  time  there  will  be  a  decrease  in 
the  double  bond  character  of  the  CN  bond.  Clearly,  the  overlap 
integral  between  the  pxr  orbitals  of  the  ring  carbon  and  of  the 
nitrogen  atom  will  show  a  strong  dependence  on  the  nitro-group 
twist  angle.  One  would  expect  that  this  overlap  integral  would 
act  as  a  gate  which  can  modulate  the  otherwise  strong 
electrophilic  properties  of  the  nitro-group.  Thus,  these  first 
order  considerations  lead  us  to  the  expectation  that  increasing 
the  nitro-group  tvist  angle  would  result  in  a  decrease  in  the  pr? 
electron  population  of  the  nitro-group. 
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Since  quadrupole  coupling  constants  are  very  sensitive  to  the 
electron  densities  local  to  the  resonant  nucleus,  one  expects  that 
the  Nitrogen-14  NQR  spectra  of  n.’.  trobenzenes  would  provide  a  good 
test  of  these  hypotheses. 

Analyses  of  the  data  reported  in  previous  ^4N  NQR  studies  of 
34  T5 

nitrobenzenes  dealing  mainly  -with  planar  molecules,  disclose 

that  of  the  three  nitrogen  occupation  numbers,  aNQ,  c?NC,  and  it, 
variations  in  the  first  reflected  variations  in  resonance  as  well 
as  inductive  effects  of  various  substituents .  (The  occupation 
numbers  were  defined  on  page  3-12.)  On  the  other  hand,  the 
occupation  number  o  ^  is  found  to  ba  insensitive  to  resonance 
effects.  Thus  of  the  Townes-Dailay  relations,  Equation  3.2-2c  is 
the  one  to  be  used  in  a  study  of  the  nitro-group  twist  since  it 
best  separates  resonance  from  inductive  effects  at  the  Nitrogen 
site. 


(aNC  **  =  (2v_  "  3  V+*  /qo 


'  (3 . 2r-2c} 

(repeated  for  convenience) 
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Figure  3. 2. 2-1  is  a  plot  of  (2v_  -  -jv+)  for  each  of  the  six 

Nitrogen-14  NQR  sites  in  monoelinic  a-TNT  vs.  the  nitro-group 

28 

twist  angle,  G ,  as  measured  by  Duke  .  Note  that  the  figure  shews 
a  very  suggestive  relation  between  the  NQR  spectroscopic 
parameters  and  the  structural  parameter, 0  .  This  relation 
inaicates  that  there  is  a  loss  of  nitrogen  pr  orbital  occupation 
of  the  order  of  1/3  of  a  milli-electron  per  degree  of  nitro-group 
twist. 
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Figure  3. 2. 2-1.  A  plot  of  the  expression  (2v  -  4/3  v  ) 

as  a  function  of  the  Nitrogen-14  lines  for  each  * 

six  sites  of  monoclinic  a -TNT  at  77°K  vs.  0,  the 
angle  of  twist  of  the  nitro-group  about  the  CN  bond 
axis. 
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The  NQR  data  used  in  Figure  3. 2. 2-1  were  obtained  at  77° K; 
the  temperature  dependence  of  the  NQR  data  is  discussed  in  Section 
3.3.  The  assignment  of  crystallographic  sites  to  NQR  sites  was 
done  by  matching  sites  with  a  higher  value  of  the  twist  angle 
with  sites  with  lower  values  of  the  occupation  number  7r .  Of 
course,  this  ensures  that  a  curve  connecting  the  points  plotted  in 
Figure  3. 2. 2-1  be  monotonic.  The  argument  that  the  correlation 
depicted  in  Figure  3.2. 2-1  is  real  finally  rests  on  the  self 
consistency  of  the  assignment  and  on  its  predictive  value.  The 
temperature  dependence  of  the  NQR  lines,  discussed  in  Section  3.3, 
will  provide  part  of  this  support.  In  the  next  section  (3.2.3)  we 
turn  to  an  NQR  study  of  other  nitrobenzenes  to  test  the  predictive 
value  of  the  relation  discovered  in  Figure  3. 2. 2-1. 

3.2.3  Structure  and  Bonding  of  TNT:  Supporting  NQR  Data 

In  the  previous  section  we  showed  how  the  details  of  the  NQR 

spectrum  of  a- TNT  could  be  accounted  for  by  considering  the  effect 

on  the  nitrogen  pir  population  due  to  the  rotation  of  the  plane  of 

the  nitro-group  about  the  CN  bond  axis.  Here  we  show  how  NQR  data 

in  other  nitrobenzenes  support  the  model  we  developed.  Table 

14 

3.2. 3-1  is  a  collection  of  N  NQR  data  for  nitrobenzenes  at  77° K. 
Each  NQR  site  is  given  a  reference  number.  Numbers  1  to  12 
correspond  to  the  twelve  sites  observed  in  the  two  phases  of 
a- TNT.  Numbers  13  to  23  are  new  data  presented  here  for  the  first 
time,  and  numbers  24  to  28  represent  all  the  nitrobenzenes  for 
which  ^N  NQR  data  was  previously  available  and  which  also  satisfy 
the  two  following  additional  requirements, 

a)  A  crystal  structure  is  available36  46 

b)  No  hydrogen  bonding  is  possible 
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TABLE  3. 2. 3-1 


NITROGEN- 14  NQR  SPECTRAL  DATA  FOR  SOME  NITROBENZENES  AT 
77°K.  ALSO  LISTED  IS  THE  ANGLE  OF  TWIST  OF  THE 
N02-GR0UP  RELATIVE  TO  THE  BENZENE  RING 


■ 

COMPOUND/SITE 

v_ 

kHz 

kHz 

e‘o0/h 

kHz 

n 

(2v_-4vi/2) 

kHz 

•i 

NITRO  GROUP 
TWIST  ANGLE 
(DEGREES) 

a  -TNT  (Monoclinic) 

1 

OAI 

792.1 

869.4 

1107.7 

0.1396 

425.0 

60.063 

2 

oAll 

801.8 

695.4 

1131.5 

0.1654 

409.7 

50.47 

3 

08 1 1 

767.7 

857.0 

1083.1 

0.1649 

392.7 

46.19 

4 

oB! 

767.3 

875.6 

1095.3 

0.1978 

367.1 

40.48 

5 

Pi 

729.7 

861.8 

1061.0 

0.2490 

310.3 

32.88 

6 

pH 

730.1 

888.1 

1073.8 

0.2929 

276.1 

22.45 

a-7NT  (orthorhombic) 

7 

oAI 

792.4 

807.8 

1106.8 

0.1362 

427.7 

59. 524 

8 

oAII 

603.2 

897.9 

1134.1 

0.1670 

409.2 

50.92 

9 

0BI1 

768.5 

858.6 

1034.7 

0.1661 

392.2 

45.70 

10 

oBt 

766.5 

875.7 

1094.8 

0.1995 

365.4 

41.32 

11 

pi 

727.4 

862.3 

1059.8 

0.2546 

305.1 

32.80 

12 

Pi  I 

732.5 

885.6 

1078.7 

0.2839 

284.2 

22.33 

U 

Trinitrooentene 

758.0 

959.4 

1141.9 

0.3518 

236.8 

14 

( TN3 ) 

752.6 

957.4 

1140.0 

0.3593 

228.7 

3-10b 

15 

734.6 

931.5 

ilio.7 

0.3545 

227.1 

16 

717.2 

907. 

1082.8 

0.3506 

225.1  J 

17 

668.6 

884. 

1035.1 

0.4162 

158.5 

S’  BENO 

18 

704.0 

832. 

1024.0 

0.2500 

298.7 

28 

19 

n-Oinitrobeniene 

811.0 

1030.0 

1227.3 

0.3569 

248.7 

10. 3C 

20 

815.7 

1037.8 

1215.7 

0.3595 

247.7 

12.4 

21 

Trinilro-«-Xylene 

759  * 

822. 

422.0 

75.2' 

22 

(TtU) 

886. 

336.7 

35.7 

23 

3.4  ,6-TrlnUro-I-t- 

752.8 

$25.1 

1051.9 

0.1375 

405.5 

-8u  tyl  «3 ,5-OMMf  Ihjrl 

Xt>  t 

Av4tllblp 

( Is33) 

24 

t 

925. 

Ul. 3 

1425.3 

0.4041 

J3J.7 

0» 

25 

P-MtrotoWefic 

911. 

119.8 

KC4.0 

QMZi 

«**  -  ? 

3.0 

26 

* 

634. 

103.7 

U47.4 

0.3259 

384.9 

11.8' 

!>•;.«  trctORtOftUrl  it  * 

874 

111.6 

1327.1 

3.3631 

2(al  .0 

JO-J  K 

j 

106.6 

1223. 7 

i  25.67 

.  / 
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FOOTNOTES  TO  TABLE  3. 2. 3-1 


a)  Reference  28 

b)  Reference  41 

c)  Reference  42 

d)  Presumed  v_  doublet 

e)  Reference  43 

f)  S.N.  Subbarao,  E.6.  Sauer  and  P.J.  Bray, 
Phys.  Letters  42A,  461  (1973) 

g)  Reference  36 

h)  Reference  34 

i )  Reference  37 

j )  Reference  38 

k)  Reference  39 

l)  Reference  35 

m)  Reference  40 


The  last  requirement  is  necessary  since  it  has  long  been  known 
that  hydrogen  bonding  in  the  solid  state  contributes  its  own 
specific  perturbations  of  NQR  spectra.  These  effects  are 
irrelevant  to  the  present  arguments  and  would  unnecessarily 
complicate  our  analysis. 


The  four  additional  compounds  whose  NQR  data  were  obtained  as 
part  of  this  study  are  the  following: 


a)  s-Trinitrobenzene  (TNB) .  Three  of  these  twelve  lines 
were  reported  by  Subbarao  and  Brady34.  six  pairs  of 

v+  and  v_  lines  are  reported  here,  in  agreement:  with  the 
multiplicity  expected  from  the  results  of  a  crystal 
structure  study*!. 

b)  m-Dinitrobenzene  (DNB) .  Subbarao  and  Brady34  using  ^5 
single  shot  spin-echo  techniques,  are  Cheng  and  Brown 
using  Nuclear  Quadrupole  Double  Resonance  techniques, 
both  report  only  one  pair  of  14n  lines.  Available 
crystal  structure  studies42,  however,  indicate  the 
presence  of  two  inequivalent  N02  groups  in  this  crystal. 
Using  the  sensitivity  enhancements  and  the  high 
resolution  capabilities  of  multiple  pulse  FT-NQR 
spectroscopy  (see  Sections  2.1  and  2.2)  we  were  able  to 
easily  resolve  the  "lines"  reported  in  the  literature^ , 35 
into  doublets  of  spacing  5.7  KHz  and  7.8  KHz, 
respectively  for  the  v+  and  v_  resonances. 

c)  Trinitro-m-xylene  (TNX) .  This  compound  was  synthesized 
for  us  by  Dr.  Prince  of  Hunter  College  of  CUNY.  It  was 
chosen  because  a  crystal  structure  study43  reports  it  to 
have  such  large  nitro-group  twist  angles,  35.7°  and  75.2°, 
that  it  should  make  a  strong  test  of  the  correlation 
introduced  in  the  previous  section. 
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(d)  2 , 4 , 6-Trinitro-l-t-Butyl-3 , 5-Dimethylbenzene  (TBDB) . 

This  analogue  of  a-TNT  has  correspondingly  bulkier 
substituents  surrounding  the  s-T.rinitrobenzene  skeleton. 
In  addition,  it  is  commercially  available  from  K&K 
Laboratories,  Plainview,  NY. 

In  addition,  to  14N-NQR  data  for  all  of  the  above  compounds 

Table  3. 2. 3-1  also  lists  the  crystallographically  determined 

nitro-group  twist  angle,  9,  when  available.  Also  listed  is  the 

4 

NQR  parameter  (2  v_  -  -j  v+)  which  according  to  Equation  3. 2- 2c 
yields  a  measure  of  the  nitrogen  pn  orbital  occupation  number. 
Figure  3. 2. 3-1  is  a  plot  of  this  NQR  parameter  vs.  the  sin  0. 

Note  that  we  use  the  sin  0  function  instead  of  simply  9  as  the 
abscissa  of  our  plot.  The  sin  0  is  a  more  appropriate  function 
here  since  we  wish  to  display  the  change  of  nitrogen  pir  orbital 
occupation,  an  effect  which  is  zero  for  9=0  and  maximum  when 
0=90°. 

Figure  3. 2. 3-1  shows  that  over  a  rather  large  range  of 

the  angle  of  twist, 9,  there  is  an  approximately  linear  variation 
in  the  double-bond  character  of  the  CN  bond,  and  that  this 
relationship  is  amply  demonstrated  by  the  NQR  data. 

3 . 3  TEMPERATURE  DEPENDENCE 

The  temperature  dependence  of  the  twenty-four  v_  and  v+ 
lines  in  a-TNT  between  77°K  and  room  temperature  are  listed  in 
Table  3.3-1  (monoclinic  phase)  and  Table  3.3-II  (orthorhombic 
phase.)  Graph  of  these  data  are  depicted  in  Figure  3.3~1  and 
Figure  3.3-2,  respectively. 
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Figure  3.2. 3-1 


4 

Plot  of  the  parameter  (2v  -  v.),  shown  in 

.  >  *T’ 

the  text  to  be  proportional  to  the  nitrogen  pu 
occupation  number  of  the  NO^  group,  vs.  the 
sine  of  the  angle  of  twist  o,  of  the  N02  group 
away  from  the  plane  of  the  benezene  ring. 
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TABLE  3.3-1 


LIST  OF  TWELVE  NITROGEN- 14  NQR  SPECTRAL  FREQUENCIES  FOR 
MONOCLINIC  TNT  AS  A  FUNCTION  OF  TEMPERATURE  FROM  ?7°K  to  321°K 


T,°K 

v,kHz 

T,°K 

v,kHz 

T,°K 

v,kHz 

T,°K 

v,kHz 

77.0 

729.7 

77.0 

730.1 

77.0 

767.3 

77.0 

767.7 

97.3 

729.4 

97.3 

729.4 

96.1 

766.5 

96.1 

766.5 

160.9 

726.0 

160.9 

726.0 

166.3 

759.2 

166.3 

759.8 

211.7 

722.1 

211.7 

722.1 

215.0 

752.4 

215.0 

754.0 

249.0 

718.9 

249.0 

718.9 

222.0 

751.4 

222.0 

753.2 

265.5 

717.0 

265.5 

717.0 

234.8 

749.6 

234.8 

751.5 

272.2 

716.3 

272.2 

716.3 

251.0 

746.9 

251.0 

749.4 

288.0 

715.1 

288.0 

715.1 

265.5 

744.3 

272.2 

746.9 

306.3 

713.4 

306.3 

713.4 

272.2 

743.0 

288.0 

744.0 

321.4 

711.6 

321.4 

711.6 

320.3 

735.1 

306.3 

741.4 

321.4 

739.2 

77.0 

792.1 

77.0 

801.8 

77.0 

857.0 

77.0 

861.8 

95.5 

790.9 

106.8 

799.8 

109.9 

855.8 

116.3 

860.5 

170.7 

778.9 

174.2 

790.2 

144.0 

852.9 

146.5 

858.3 

217.7 

769.5 

221.0 

782.4 

177.4 

849.9 

180.6 

855.7 

252.3 

761.8 

253.6 

776.5 

195.0 

348.2 

198.2 

854.2 

265.5 

758. 3 

265.8 

774.1 

224.5 

845.3 

225.8 

851.7 

272.2 

756.7 

272.2 

772.9 

237.4 

844.1 

240.4 

850.3 

288.0 

753.2 

288.0 

770.0 

265.8 

840.6 

265.8 

847.6 

306.3 

748.6 

306.3 

766.5 

272.2 

839.8 

272,2 

846.8 

321.6 

744.6 

321.6 

763.5 

288.0 

838.1 

288.0 

845.1 

306.3 

835.8 

306.3 

842.6 

321.6 

833.7 

320.6 

841.0 

77.0 

869.4 

77.0 

875.6 

77.0 

888.1 

77.0 

895.4 

119.4 

867.0 

124.9 

872.0 

130.0 

885.8 

132.4 

889.7 

148.9 

864.0 

151.1 

868.9 

155.2 

884.0 

157.2 

885.5 

132.6 

860.9 

184.5 

864.6 

188.0 

881.3 

188.0 

880.6 

201.9 

857.6 

204.9 

862.0 

209.9 

879.3 

209.9 

876.9 

228.5 

853.8 

229.5 

858.4 

233.8 

877.0 

232.8 

872.4 

240.4 

852.2 

242.7 

856.6 

248.3 

875.6 

246.0 

870.1 

265.5 

847.5 

265.5 

852.9 

266.3 

873.4 

256.0 

365.5 

272.2 

846.8 

272.2 

851.9 

272.2 

872.6 

272.2 

364.2 

288.0 

844.1 

288.0 

849.6 

288.0 

871.1 

288.0 

360.9 

306.3 

840.7 

306.1 

347.6 

306.1 

865.2 

306.1 

856.6 

321.6 

837.5 

321.6 

844.0 

321.1 

867.3 

321.1 

853.3 
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TABLE  3. 3- 1 1 


LIST  OF  TWELVE  NITROGEN-14  NQR  SPECTRAL 
FREQUENCIES  FOR  ORTHORHOMBIC  TNT  AS  A 
FUNCTION  OF  TEMPERATURE  FROM  77° K  TO  279.6° K 


T,°K 

v,kHz 

- 1 - 

T,°K 

v,kHz 

T,°K 

v,kHz 

T,°K 

v,kHz 

77.0 

eh 

77.0 

732.5 

77.0 

766.5 

77.0 

768.5 

81.8 

El 

81.8 

732.4 

83.2 

766.5 

83.2 

768.5 

87.4 

KpH 

88.8 

732.4 

90.9 

766.2 

90.9 

768.1 

106.8 

EHfl 

106.8 

731.8 

116.3 

763.8 

116.3 

765.6 

127.4 

725.3 

113.2 

731.4 

126.6 

762.8 

126.6 

764.8 

156.7 

723.4 

129.3 

730.3 

144.8 

761.0 

144.8 

762.4 

170.0 

722.5 

141.4 

729.6 

163.1 

753.8 

163.1 

760.1 

187.8 

721.4 

173.4 

727.3 

175.6 

757.3 

175.6 

758.5 

205.2 

720.0 

189.5 

726.1 

179.6 

756.9 

179.6 

758.0 

214.0 

719.3 

206.9 

724.8 

192.0 

755.3 

192.0 

756.3 

225.3 

718.4 

215.7 

724.0 

208.7 

753.4 

208.7 

754.1 

233.9 

717.1 

226.5 

723.1 

217.2 

752.1 

217.2 

752.9 

279.6 

713.3 

236.9 

721.9 

227.8 

750.7 

227.8 

751.5 

279.6 

718.0 

240.9 

749.4 

240.9 

749.4 

279.6 

743.7 

279.6 

743.7 

77.0 

792.4 

77.0 

803.2 

77.0 

858.6 

77.0 

862.3 

84.1 

792.3 

85.7 

803.1 

81.3 

858.3 

81.3 

862.1 

96.1 

791.2 

101.8 

801.4 

83.2 

858.6 

83.2 

862.2 

119.7 

787.3 

123.7 

798.4 

90.9 

853.0 

90.5 

861.9 

150.6 

782.6 

153.5 

794.0 

104,9 

857.3 

108.7 

861.1 

165.3 

780.0 

165.3 

792.6 

128.8 

855.2 

131.9 

860.0 

182.8 

776.9 

185.0 

789.0 

131.7 

854.9 

138.2 

859.5 

197.0 

774.0 

199.2 

786.8 

153.0 

852.9 

1:35.5 

858.4 

210.7 

771.3 

212.5 

784.4 

163.3 

852.0 

167.7 

857.7 

220.2 

769.1 

221.7 

782.6 

172.4 

851.1 

175.4 

857.1 

229.0 

767.3 

229.0 

781.0 

183.8 

849.7 

179.3 

856.7 

242.7 

764.7 

244.2 

778.6 

191.0 

849.0 

187.5 

856. 2 

279.6 

755.6 

279.6 

771.8 

193.5 

348.8 

196.0 

855.? 

204.2 

847.7 

208.9 

854.6 

217.7 

846.0 

221.0 

853.5 

219.0 

846.2 

235.3 

852.3 

220.2 

848.0 

231.6 

844.7 

350.5 

233.1 

844.6 

279.6 

848.0 

246,2 

842.9 

256.1 

842.0 

279.6 

839.5 
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TABLE  3. 3-1 1  (CONCLUDED) 

LIST  OF  TWELVE  NITROGEN-14  NQR  SPECTRAL 
FREQUENCIES  FOR  ORTHORHOMBIC  TNT  AS  A 
FUNCTION  OF  TEMPERATURE  FROM  77°K  TO  279.6° K 


T,°K 

v,kHz 

T,°K 

v.kHz 

T,°K 

v.kHz 

T,°K 

v,kHz 

77.0 

867.8 

77.0 

875.7 

77.0 

885.6 

77.0 

897.9 

80.4 

867.8 

80.1 

875.6 

80.4 

885.5 

80.6 

897.8 

84.3 

867.8 

85.5 

875.6 

86.5 

885.5 

87.6 

897.5 

91.9 

867.5 

92.8 

875.1 

98.7 

884.9 

119.0 

892.5 

134.6 

864.0 

95.2 

875.0 

103.2 

884.8 

142.3 

890.0 

141.4 

863.0 

137.2 

870.5 

115.6 

884.0 

155.2 

888.0 

157.9 

861.0 

143.8 

869.6 

146.2 

881.8 

165.3 

886.4 

170.0 

860.0 

160.6 

867.6 

163.1 

880.5 

177.9 

884.7 

181.8 

858.4 

177.6 

865.6 

175.1 

879.5 

184.5 

883.1 

189.7 

857.6 

192.2 

863.3 

182.6 

879.0 

198.7 

880.6 

207.7 

855.1 

194.0 

863.2 

196.5 

877.8 

205.7 

879.4 

221.5 

853.6 

200.0 

862.3 

201.5 

877.6 

237.6 

873.4 

248.5 

849.5 

201.5 

861.9 

214.2 

876.4 

253.8 

870.2 

258.4 

848.1 

210.0 

860.5 

224.5 

875.6 

279.6 

864.7 

279.6 

844.6 

222.7 

859.1 

226.0 

857.3 

224.5 

858.9 

227.3 

875.3 

229.3 

858.3 

227.5 

875.1 

236.4 

857.1 

237.6 

874.2 

250.3 

855.1 

239.4 

874.0 

262.5 

353.4 

253.8 

872.9 

279.6 

850.5 

279.6 

870.0 

x  800 
o 

x: 
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Figure  3.3-1.  Temperature  Dependence  of  the 
Nitrogen-14  NQR  Spectrum  of 
MONOCLINIC  TNT 
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Thermal  motions  present  in  all  crystals  occur  at  frequencies 
that  are  much  greater  than  NQR  transitions,  thus  the  electric 
field  gradients  actually  measured  by  NQR  techniques  are  time- 
averaged  over  these  motions.  Larger  librational  amplitudes  are 
generally  more  efficient  at  this  averaging  process,  so  that, 
typically,  NQR  frequencies  are  monotonically  decreasing  functions 
of  the  temperature.  Figures  3.3-1  and  3.3-2  clearly  demonstrate 
this  feature.  It  is  generally  a  good  assumption  to  describe 
the  torsional  motions  by  an  Einstein  model,  and  then  the  root 
mean  square  of  the  amplitudes  about  the  x,y,  and  z  inertial 
axes  of  the  nitro-greup  are  given  by 


<A0 . 2>  =  — ~ 
i  ,-_2 


4it  I .  v  . 

l  l 


.1  +  _ I _ 

2  exo(hv./kT)  -1 


;  i-x,y, z 


(3.3-1) 


« 

where  1^  is  the  effect-)  ve  moment  of  inertia  about  the  i 
principal  axis  and  is  che  torsional  frequency.  To  the 
extent  that  the  principal  axes  of  the  inertial  tensor  coincide 
with  the  electric  field  gradient  tensor  principal  axes,  the 

change  in  the  coupling  constant  and  asymmetry  parameter  can 

.  ,  44,45 

be  expressed  as 


(e2qQ) 


_  i 
1  2 


and 


(e  'qQ)T=0 


(T)  (e  qQ), 


n,T=0)  (e2qQ).t=0 


3 


<A0  2>  +  <A0  2> 
x  y 


+  f 

f  2  2 

<A0  >  -  <A0  > 

-2  <A0  2> 

2'.\ 

y  x 

z 

.  j 

■ 

’ 

r 

<A3  2>  +  <A0  2> 

+  ^ 

<A0  2> 

-  <A0  2> 

x  y 

2 

y 

X 

.  A 

( 

'i 

j 

(3.3-2 

(3.3-3) 
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An  analysis  of  the  temperature  dependence  of  both  phases  of 

a -TNT  shows  that  the  main  features  of  the  data  can  be  explained 

bv  considering  only  the  torsional  motion  of  each  N0_  group  about 

/  2 

the  C-N  bond,  i.e.  by  considering  only  the  effect  of  <A0  >  and 

2  2  x 
assuming  that  <A0  >  =  <A0  >  -  zero.  Here,  we  are  labelling  the 

y  z 

axes  according  to  the  convention  given  in  Figure  3. 2. 1-1. 

With  these  assumptions,  Equation  3.3-2  can  be  rewritten  as, 


<A0 

x 


2 


> 


<A9> 


2 

RMS 


(e2qQ)T=0  - (e2qQ)T 
2hv+ 


(3.3-4) 


Figures  3.3-1  and  3.3-2  show  that  at  77°K  there  is  little 

temperature  dependence,  so  we  make  only  a  small  error  by  taking 

the  "T=0"  values  at  this  temperature.  Table  3.3-III  shows  the 

components  of  the  electric  field  gradient  (EFG)  tensor  q  ,  q  ,  and 

yy 

q  at  both  77°K  and  321°K  for  monoclinic  TNT.  The  EFG  tensor 
zz 

components  were  computed  from  the  frequencies  using  the  relations. 


! 

X 

ro|<N 

II 

h  (v+  +  v_) / (a2qoQ) , 

(3.3-5) 

\ 

q 

__  ^ 

q  (1  +  n) /2, 

(3.3-6) 

} 

yy 

^ZZ 

f 

q 

nxx 

=  - 

qzz  (1  -  n)/2. 

(3.3-7) 

, .  In  the  fourth  column  of  Table  3.3-III  we  compute  Aq,  the 

I  change  in  EFG  component  due  to  the  rise  in  temperature.  Note  that 

if  an  average  of  Aq  is  taken  over  molecules  I  and  II,  as  is  done 

|  in  the  fifth  column  of  the  Table,  it  becomes  clear  that  in  fact 

the  major  change  is  due  solely  to  changes  in  q  .As  proof  for  this 

z  z 

|  assertion,  note  that  the  relations  (Aq  |  *  j  Aq  |  *  j  %  q  {  hold 

f  xx  y  y  &  ** 

reasonably  well,  as  can  be  ascertained  by  comparing  tne  entries 
[  joined  by  the  double  pointed  arrow. 
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TABLE  3.3-III 


THE  COMPONENTS  qvv,  q  .  and  q77  OF  THE  ELECTRIC 

aA  j fy  Zi- 

FIELD  GRADIENT  TENSOR  OF  MONOCLINIC  TNT  AT  77°K 
AND  AT  321° K. 


- - - 

Aqj+Aqn 

2 

TNT  SITE 

T=77°K 

T=321 °K 

Aq 

-q 

^xx 

-q 

Mxx 

oAI 

oAII 

476.1 

472.2 

434.4 

449.1 

42.1  \ 

23.1  f 

32.6  «->• 

26.6 

oBII 

oBI 

452.3 

439.3 

425.8 

417.4 

26.5  i 
21.9  / 

24.2 

19.3 

Pi 

pH 

398.4 

381.4 

388.1  ■ 
370.6 

10.3  \ 
10.8  } 

10.5 

13.1 

"qyy 

-q 

-JDL 

oAI 

oAII 

631.2 

659.3 

620.3 

628.8 

10.9  ) 
30.5  J 

-  20.7  ^ 

26.6 

oBII 

oBI 

630.9 

656.0 

622.8 

635.3 

8.1  | 
20.7  ( 

14.4  -M- 

19.3 

Pi 

pH 

662.6 

697.4 

647.0 

628.0 

15.6 

15.4 

►  15.5  ++ 

13.1 

qzz 

oAI 

1107.7 

1054.7 

53.0 

\  53.3 

oAII 

1131.5 

1077.9 

53.6 

oBII 

1083.1 

1048.6 

34.5 

|  38.6 

oBI 

1095.3 

1052.7 

42.6 

Pi 

1061.0 

1035.1 

25.9 

|  26.1 

pH 

1078.8 

1052.6 

26.2 

/ 


Table  3.3-IV  is  analogous  to  Table  3.3-III  in  that  it 
lists  the  EFG  components  for  the  orthorhombic  phase  of  TNT 
for  both  77°K  and  279. 6°K.  Here  also  it  is  approximately 
true  that  |  Aqxx  |  *|  A<g^.yj  =»  |-|  Aqzz,|  validating  our  assumption 
that  A9  *  A0  *  0. 

y  z 

Finally,  in  Table  3.3-V  the  root  mean  square  angles  of 
nitro*-group  libration  ,  <AQX>  RMS,  are  computed  from  Equation  3.3-4 
and  listed  together  with  R.C. J.  Duke's  values  of  the  equilibrium 
angle  of  nitro-group  twist.  Note  that  the  angular  displacements 
are  considerable  in  magnitude  in  all  cases,  and  are  larger  for 
larger  values  of  0.  This  can  be  understood  as  an  indication 
that  cases  with  a  larger  value  of  9,  and  therefore  with  a 
weaker  double-bond  character,  have  a  correspondingly  less  stiff 
torsional  spring  constant.  Thus  larger  torsional  angles  are 
reached  at  any  given  temperature.  The  implications  of  these 
considerations  to  the  thermal  break  up  of  the  a-TNT  molecule 
deserve  further  study. 
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TABLE  3. 3- IV 

THE  COMPONENTS  q  .  q  ,  and  q  OF  THE  ELECTRIC 
xx  yy  zz 

FIELD  GRADIENT  TENSOR  OF  ORTHORHOMBIC  TNT  AT  77°K 
AND  AT  279. 6°K. 


Aqj+Aqjj 

2 

K 

77°K 

279. 6°K 

Aq 

_qxx 

-q 

Mxx 

oAI 

478.0 

444.4 

33.6  1 

26.7  -w 

20.8 

oAH 

472.3 

452.6 

19.7  f 

oBII 

452.3 

431.9 

20.4  \ 

17.0  *■+ 

15.3 

oBI 

438.2 

424.6 

13.6  / 

Pi 

PH 

395.0 

386.3 

386.4 

377.3 

8.6  1 
9.0  / 

8.8 

9.7 

-Q 

Hyy 

-q 

JH. 

oAI 

628.8 

622.4 

6.4 

\  14.9  ** 

20.8 

OAII 

661.7 

638.4 

23.3 

oBII 

632.5 

623.5 

9.0 

}  13.7 

15.3 

OBI 

656.6 

638.2 

18.4 

PI 

PII 

664.8 

692.5 

654.8 

681 .3 

10.0 

11.2 

|  10.6  «-*■ 

9.7 

qzz 

qzz 

oAI 

1106.8 

1066.8 

43.0 

|  41.6 

oA  II 

1134.1 

1091.0 

43.1 

oBII 

1084.7 

1055.5 

29.2 

f  30.6 

oBI 

1094.8 

1062.8 

32.0 

Pi 

1059.8 

1041.2 

18.6 

\  19.3 

PH 

1078.7 

1058.7 

20.0 
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TABLE  3.3-V 

THE  ROOT  MEAN  SQUARE  LIBRATIONAL 
ANGLE  <A6x>rms  AT  ROOM  TEMPERATURE 
FOR  MONOCLINIC  (44°C)  AND  ORTHORHOMBIC 
(2.6°C)  TNT.  ALSO  SHOWN  IS  THE 
CRYSTALLOGRAPHIC  NITRO  GROUP  TWIST 
ANGLE  DEFINED  IN  THE  TEXT. 


SITE 

<4Vrms 

0a 

MONOCLINIC  TNT 

oAI 

10.2° 

60.06° 

oAII 

10.2° 

50.47° 

oBII 

8.2 

46.19° 

oBI 

9.1° 

40.48° 

Pi 

7.1° 

32.88° 

pH 

• 

o 

o 

22.45° 

ORTHORHOMBIC  TNT 

oAI 

8.8° 

59.52° 

oAII 

9.0° 

50.92° 

oBII 

7.6° 

45.70° 

oBI 

7.9° 

41.32° 

Pi 

6.0° 

32.80° 

PII 

6.1° 

22.33° 

a)  From  R.C.U.  Duke,  Private  Communication  (1982) 
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4.0  NQR  STUDY  OF  HMX  AND  RDX46 
4 . 1  INTRODUCTION 


1,3, 5,7-  Tetranitro-l,3,5,7-Tetraazacyclooctane  (HMX),  shown 

in  Figure  4.1-1,  has  been  shown  to  exist  in  several  crystal 

phases.  For  example,  T.B.  Brill  and  coworkers  have  used  laser 

Raman  spectroscopy  to  diagnose  polymorphic  interconversions  as  a 

47  48  4Q 

function  of  temperature  and  pressure  '  '  ' .  The  stable  form  at 

room  temperature  is  labelled  8 -HMX,  while  <$ -HMX  is  the 

thermodynamically  favorable  form  under  the  conditions  of  high 

49 

pressure,  for  example  those  in  a  rocket  motor 

Questions  remain  about  the  rate  of  the  p->6  solid-solid 

transformation  and  mechanism  of  the  solid  phase  events  up  to  and 

including  the  phase  transition.  These  problems  relate  directly  to 

modelling  of  and  performance  in  the  combustion  environment.  While 

Raman  spectroscopy  is  useful  for  distinguishing  the  HMX 
47 

polymorphs  and  for  specifying  some  of  the  events  leading  to 

50 

degradation  m  inorganic  materials,  it  has  not  proven  to  be  as 

useful  in  elucidating  the  detailed  events  leading  to  degradation 

in  HMX.  This  difficulty  appears  to  result  from  the  ready 

dissipation  of  thermal  energy  among  the  numerous  internal 

stretching  and  bending  modes  of  the  HMX  molecule.  The  vibrational 

spectrum  of  B-HMX  shows  gradual  non-specific  changes  with 

47 

temperature  up  to  the  phase  transition. 
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Figure  4.1-1.  The  HMX  Molecule 


i 
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Information  about  molecular  d -ramies  in  KMX  should  be 


obtainable  if  a  spectroscopic  technique  is  chosen  which  samples  on 

a  time  scale  of  molecular  rather  than  atomic  motion.  The 

mechanisms  of  solid  phase  events  should  be  definable  by  this 

approach.  Nitrogen-14  nuclear  quadrupole  resonance  spectroscopy 

is  just  such  a  method,  and  it  was  employed  for  these  studies.  NQR 

frequencies  aid  in  understanding  the  electron  distribution  in 
2 

molecules,  while  the  temperature  dependence  of  the  resonance 

frequencies  has  been  used  to  cast  light  on  molecular  motion  in 

...  51,52,53,54 

organic  solids. 

14 

We  report  here  an  analysis  of  the  N  NQR  data  for  the  amine 

nitrogen  atoms  of  HMX  in  terms  of  the  electronic  structure  and 

iibrational  motions  of  the  molecule.  The  results  are  interpreted 

in  light  of  the  intermolecular  cohesive  forces  that  exist  in  the 

14 

crystal  lattice  of  B-HMX.  The  N  NQR  data  for  the  amine  nitrogen 
atoms  in  RDX  (hexahydro-1 . 3 , 5-trinitro-s-triazine)  ,  which  is  the 
six  membered  analog  of  HMX,  are  also  reported  here  for  purposes  of 
comparison. 

4 . 2  EXPERIMENTAL 

Samples  of  HMX  taken  from  the  same  batch  were  used  without 
further  purification  in  these  studies.  The  particle  sise  was 
175 ym . 


4-3 


BLOCK 


The  nuclear  quadrupole  resonance  measurements  were  made  using 

a  pulsed  spectrometer  whose  details  have  been  described 
55 

elsewhere.  Heat-resistant ,  high-Q,  NQR  coils  were  constructed  by 
wrapping  Cu  wire  around  a  cylindrical  sample  bottle  and  coating 
the  coil  with  high  temperature  resin  to  ensure  that  the  coil  would 
retain  its  geometric  integrity  throughout  the  temperature  range. 

Measurements  of  the  frequency  as  a  function  of  temperature 

were  obtained  by  using  a  strong  off-resonance  comb  of  RP  pulses 

4 

(SORC) ,  2-6  KHz  off-resonance.  The  frequency  was  determined  by 

directly  measuring  the  period  of  oscillation  from  the  transient 

4 

NQR  signal  produced  between  pulses  after  10  coadds  on  a  PAR  4202 
signal  averager.  This  frequency  measurement  requires  only  about 
10  sec.  and  produces  an  error  of  +lKHz. 

Temperature  variations  below  300° K  were  achieved  by  allowing 
very  slow  heat  loss  to  occur  from  the  sample.  The  NQR  coil  was 
mounted  in  a  thick-walled  Cu  can.  The  can  was  placed  in  a  bath  of 
liquid  N^.  The  temperature  of  the  sample  was  measured  with  a 
Fluke  digital  thermometer  and  an  Omega  stainless  steel  sheathed 
thermocouple.  The  thermocouple  was  wrapped  with  Teflon  tape  and 
placed  directly  into  the  HMX  sample.  The  Teflon  tape  prevented 
any  possible  reaction  between  the  sample  and  the  thermocouple. 

After  the  temperature  equilibrated  at  77°  K,  it  was  removed  from 
the  bath,  placed  in  a  Dewar  flask,  and  tightly  insulated.  The 
frequency  measurements  were  made  every  10°K  by  allowing  the  sample 
to  warm  slowly  via  thermal  loss.  The  temperature  of  the  sample 
was  monitored  continuously  and  all  mesurements  were  made  while  the 
sample  was  within  0.5°K  of  the  desired  temperature. 
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For  the  frequency  measurements  above  room  temperature,  a  2  liter 
Dewar  flask  was  filled  with  silicone  oil.  An  Omega  model  402 
temperature  controller  coupled  to  an  immersion  heater  and  an 
iron-cons tantan  thermocouple  was  used  for  temperature  regulation. 

A  mechanical  stirrer  prevented  temperature  gradients  in  the  bath. 

The  temperature  of  the  sample  was  monitored  directly  by  the  method 
described  above.  A  vial  of  HMX  was  placed  in  the  NQR  coil  and 
immersed  in  the  silicone  oil.  The  sample  was  heated  in  a  stepwise 
fashion  through  the  temperature  range  of  310°K-4  30°K.  Frequency 
measurements  were  taken  in  approximately  10°  K  increments  by 
heating  at  a  rate  of  about  one  degree  per  minute  to  the  desired 
temperature  and  then  equilibrating  for  5  nun. 

The  variable  temperature  experiments  were  repeated  several 
times,  and  a  fresh  sample  of  HMX  was  used  each  time.  However,  the 
resonance  frequencies  were  observed  to  be  the  same  whether  the 
temperature  run  began  at  310° K  or  375° K.  As  long  as  the 
temperature  was  kept  below  the  solid  phase  transformation  at 
4  36°K,  the  resonance  frequencies  upon  cooling  to  a  specific 
temperature  were  the  same  as  those  obtained  upon  heating  to  that 
temperature. 

4.3  RESULTS  DISCUSSION 

nuclei  (1=1)  yields  three  NQR  signals  when  •  In  order 
to  assign  the  transitions  correctly,  was  measured  at  77°K  and 
298  °K  in  HMX.  The  NQR  data  for  a  8 -HMX  over  the  temperature  range 
of  77-426 °K  are  compiled  in  Table  4.3-1.  These  signals  arise  from 
the  amine  nitrogen  atoms  of  the  ring.  No  attempt  was  made  to 
obtain  data  for  the  iruch  1  awer- frequency  nitro-group  nitrogen 

14  9  -4 

sites.  N  NQR  data  at  77  K  for  the  amine  nitrogen  atoms  in  RDX 
are  given  in  Table  4.3-II. 
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TABLE  4.3-1 


14N  NQR  DATA  FOR  THE  AMINE  NITROGEN  ATOMS  OF  b-hmx 


T,  °K 

v» 

KHz 

n 

e^Oq,  kHz 

v+ 

v  _ 

77a 

5336 

3745 

0.5256 

6054 

5091 

3628 

0.5030 

5813 

83 

5336 

3745 

0.5256 

6054 

5090 

3628 

0.5031 

5812 

93 

5334 

3744 

0.5254 

6053 

5089 

3628 

0.5028 

5811 

103 

5333 

3744 

0.5252 

6051 

5088 

3628 

0.5025 

5811 

113 

5331 

3744 

0.5246 

6050 

5087 

3628 

0.5022 

5810 

123 

5330 

3744 

0.5244 

6049 

5086 

3628 

0.5020 

5809 

135 

5328 

3744 

0.5238 

6048 

5085 

3628 

0.5017 

5809 

143 

5328 

3744 

0.5238 

6048 

5085 

3628 

0.5017 

5809 

153 

5326 

3743 

0.5237 

6046 

5084 

3628 

0.5014 

5808 

163 

5325 

3743 

0.5234 

6045 

5083 

3628 

0.5011 

5807 

173 

5324 

3742 

0.5235 

6044 

5082 

3627 

0.5012 

5806 

183 

5322 

3742 

0.5229 

6043 

5080 

3627 

0.5006 

5085 

193 

5320 

3742 

0.5224 

6041 

5079 

3627 

0.5003 

5805 

203 

5318 

3742 

0.5219 

6040 

5078 

3627 

0.5001 

5803 

213 

5317 

3742 

0.5216 

6039 

5077 

3627 

0.4998 

5802 

223 

5315 

3/42 

0.5210 

6038 

5075 

3627 

0.4992 

5801 
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TABLE  4.3-1  (CONTINUED) 


TABLE  4.3-1  (CONCLUDED) 


14N  NQR  DATA  FOR  THE  AMINE  NITROGEN  ATOMS  OF  B-HMX 


T,  °K 

KHz 

n 

e2Qq,  kHz 

V+ 

v_ 

405 

5273 

3730 

0.5142 

6002 

5040 

3611 

0.4955 

5767 

415 

5272 

3728 

0.5147 

6000 

5038 

3610 

0.4954 

5765 

426 

5270 

3726 

0.5149 

5997 

5032 

3607 

0.4948 

5759 

a)  observed  at  1591,  1462  KHz 

b)  observed  at  1564,  1441  KHz 


2 

The  nuclear  quadrupole  coupling  constant,  e  Qq/h,  and  the 
asymmetry  parameter,  n  /  can  be  calculated  from  the  frequencies  of 
v+  and  v_  by  using  Eqs.4.3-1  and  2.^ 

e2Qq/h  =  2/3  (v_  +  v+)  (4.3-1) 

(v+  -  v  ) 

n  =  3  -  (4.3-2) 

(v+  + 

The  coupling  constant  measures  the  largest  component  of  the 
electric  field  gradient  (EFG)  tensor,  while  q  measures  the 
departure  of  the  EFG  from  axial  symmetry  about  this  component. 

The  four  amine  nitrogen  atoms  reduce  to  two  pairs  because  of 
the  Ch  symmetry  in  the  8-HMX  molecule  shown  in  Figure  4.3-1. 

These  two  pairs  (N2  2,  and  N3  3,)  produce  NQR  signals  whose 
coupling  constants  differ  by  241  kHz  at  77° K.  This  difference  is 
about  four  percent  of  the  overall  EFG  and  is  larger  than  that 
usually  observed  for  simple  crystallographic  inequivalence  of 
sites  in  organic  molecular  crystal.  However,  the  difference  is 
smaller  than  usually  arises  from  chemical  inequivalence.  For 
convenience  the  amine  atoms,  3,  of  8 -HMX  are  referred  to  here 
as  axial  amines,  while  N2  2,  are  labeled  the  equatorial  amine 
atoms  in  reference  to  the  positions  of  the  N02  groups  bound  to 
these  nitrogens  with  respect  to  the  plane  of  the  ring. 
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Electronic  Effects;  Assignment  of  the  nuclear  quadrupole 
coupling  constants  to  the  axial  and  equatorial  sets  of  the 

nitrogen  atoms  can  be  made  by  employing  the  Townes-Dailey 

39  S7 

formalism  and  the  crystal  structure  of  g-HMX.  .  In  the 

Townes-Dailey  model  the  diagonal  components  of  the  EFG  tensor  are 

described  in  terms  of  the  electron  populations  of  the  atomic 

orbitals  on  the  nitrogen  atom. 

The  molecular  structure  of  HMX  reveals  that  the  C-NNO„  unit 

57  2  2 

of  the  molecule  is  nearly  planar.  As  a  result,  the  a -bonding 

orbitals  on  the  amine  nitrogen  atoms  can  be  constructed  assuming 

sp  hybridization.  It  is  helpful  to  consider  the  orientation  of 

the  EFG  principal  axes  when  constructing  these  hybrids  in  order 

that  the  convention  of  |q  \ > | q  I  >  |  a  I  is  followed.  The  proper 

z  z  y  y  XX 

coordinate  system  usually  can  be  chosen  on  the  basis  of  the 

relative  electronegativities  of  the  atoms  or  an  MO  calculation,  if 

available.  Unambiguous  assignment  of  the  principal  axes  requires 

a  single  crystal  Zeeman  study  or  other  experiments.58  However,  in 

the  HMX  case  the  orientation  of  the  principal  axes  can  be 

determined  with  some  confidence  without  such  studies.  This  is  due 

2 

to  the  fact  that  the  sp  hybridized  amine  nitrogen  atoms  contain  a 
pair  of  elctrons  having  ir-bonding  symmetry  perpendicular  to  the 
C2NN02  plane.  Placing  the  Z  axis  of  the  EFG  coincident  with  this 
electron  pair  is  consistent  with  other  studies  involving 
substituted  nitrogen  atoms  having  sp  hybridization,  -s  urea, 59 

thiourea,6^  N-methy  1-pyrrole, ^  and  hydroxypyrimidines.0,L 
The  Y  axis  coincides  with  the  N-N  bond  and  the  X  axis  is 
perpendicular  to  the  YZ  plane.  Figure  4„3-2  pictorializes  this 
axis  system  which  was  also  used  to  construct  the  hybrid  atomic 
orbitals  listed  in  Table  4.3-III. 
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Figure  4.3-2.  The  approximate  orientation  of  the  principal 
ax"  of  the  EFG  tensor  of  the  amine  atoms. 

Z  is  perpendicular  to  the  molecular  plane. 


TABLE  4.3-III 


THE  HYBRID  ORBITALS*  OF  THE  AMINE  NITROGEN  ATOMS 
ACCORDING  TO  THE  AXiS  SYSTEM  GIVEN  IN  FIGURE 
THE  ORBITAL  OCCUPATION  NUMBERS  ARE  ALSO  GIVEN 


e  -  cot  (y)  *Ntt 
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The  Townes-Dailey  Eq.  4.3-3  relates  the  diagonal  components 
of  the  EFG  tensor,  q^,  to  the  2p  orbital  populations  on  nitrogen, 
N  . 


e2Qqii/h  =  e2Qqp/h 


N  .  N 

Np  _  -Zl  ■  Pk 

Fi 


x,y,z  =  i,j,k  (4.3-3) 


2 

e  Qqp/h  is  the  coupling  constant  produced  by  a  single  2p  electron 
on  nitrogen.  When  the  wavefunctions  in  Table  4.3-III  are 
substituted  into  eq  (4.3-3),  one  obtains  Eq.4.3-4  and  4.3-5  which 
describe  the  coupling  constant  and  asymmetry  parameter,  n ,  for  the 
amine  nitrogens  in  terms  of  the  hybrid  orbital  populations,  a,b, 
and  c.  * 


’  \ 

e2Qq/h 

n  3 

e  Qqp/h 

k  4 

2 

'NC 


> 

* 

.  2 
l-cot  v 

. 

r  ' 

e2Qq/h 

e2QqD/h 

. 

aNu  ”  Sc 


(4.3-4) 


(4.3-5) 


Eq  4.3-4  contains  the  dependence  of  the  EFG  parameters  on  y ,  which 
is  one  half  the  CNC  angle  in  HMX.  Inclusion  of  this  refinement  in 
the  Townes-Dailey  model  hardly  affects  the  results  in  HMX  because 
Y°  61.9°  at  the  axial  amine  and  61.2°  at  the  equatorial  amine. ^ 
Thus,  eliminating  the  angle  term  reduces  Eq  4.3-4  to  Eq  4.3-6. 


4-15 


BLOCK 


e2Qq/h 

e2Qqp/h 


n  =  c. 


NC 


-  b. 


NN 


(4.3-6) 


In  Eq.  4.3-6  it  is  assumed  that  the  population  of  the  orbital 

involved  in  the  NN  bond  is  less  than  that  in  the  NC  bond.  This 

difference  is  both  chemically  reasonable  and  consistent  with  ab 

6  2 

initio  MO  calculations  on  HMX. 


The  atomic  coupling  constant  for  a  single  p  electron  on  a 
2 

nitrogen  atom,  e  Qq  /h,  has  been  taken  to  be  8.4  MHz  although  the 

Jr 

exact  value  cannot  be  determined  experimentally.  Uncertainty  in 
this  value  is  of  importance  when  absolute  orbital  populations  are 
sought,  but  is  not  vital  when  making  relative  comparisons  as  will 
be  done  here.  Nevertheless,  Eqs.  (4.3-5)  and  (4.3-6)  still 
contain  three  unknowns,  a^ ,  b^ ,  and  cNC  but  only  two 
experimental  observables,  e2Qq/h  and  n.  An  assumption  about  one 
of  the  unknowns  must  be  made.  If  a  value  of  l.lOe  is  chosen  for 
cM„  based  on  previous  NQR  analyses,2^  then  an  indication  of  the  a 
and  tr  character  in  the  N-N  bond  can  be  deduced.  The  resulting 
values  of  aNlI  and  bNN  are  compiled  ;.n  Table  4.3-IV.  The  higher 
coupling  constant  in  HMX  can  be  seen  to  be  associated  with  a 
slightly  higher  population  for  the  amine  p2  orbital  involved  in 
the  ttNN  bond,  (i.e.,  less  it  character)  and  slightly  lower 
population  in  the  bond,  while  the  reverse  is  true  of  the  lower 
coupling  constant.  The  neutron  diffraction  data  for  8- HMX  raveal 
that  the  N-N  bond  is  longer  for  the  axial  amine?  (1.373(5)  A] 
than  the  equatorial  amines  (1.354(5)  A],  Therefore,  the  larger 
coupling  constant  arises  from  the  axial  amines  and  the  smaller 
value  is  due  to  the  equatorial  amines. 


57 
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TABLE  4.3-IV 


THE  HYBRID  ORBITAL  OCCUPATIONS  OF  THE  AMINE 
NITROGEN  ATOMS  ASSUMING  A  POPULATION  FOR  THE 
N-C  BOND  OF  1.10  ELECTRONS  AND  THE  NQR  DATA 
AT  77° K 


aNrr 

bNN 

CNC 

HMX 

N2  2  "(equatorial ) 

1,675 

0.76 

(1.10) 

N3  3  "(axial) 

1.69 

0.73 

(1.10) 

RDX 

(average) 

1.65 

0.67 

(1.10) 
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The  difference  in  the  axial  and  equatorial  amine  coupling 
constants  may  originate  in  the  intermolecular  interactions  of  the 
6-HMX  lattice.  An  examination  of  the  difference  between  the  amine 
sites  reveals  that  the  axial  fragment  is  engaged  in  a  close 

interaction  with  an  oxygen  atom,  0V  of  an  axial  NO-  group  in  the 
neighboring  molecule  (Figure  4.3-3).  '  The  carbon  atoms  and 

C2  are  electrostatically  attracted  to  the  neighboring  0^  atom 
while  the  amine  nitrogen  N3  is  repelled  by  it.  This  repulsion 
occurs  in  the  electron  region  of  and  could  be  responsible  for 
the  slight  lengthening  of  the  axial  NN  bond  relative  to  the 
equatorial  NN  bond,  where  no  such  intermolecular  interactions 
takes  place. 

The  six  membered  ring  analog  of  HMX ,  known  as  RDX,  has  the 
stoichiometry  C3HgN3 (N02) 3>  The  crystal  contains  three 
crystallographically  inequivalent  nitrogen  atoms  64  which  give 
rise  to  the  NQR  data  in  Table  4.3-II.  The  average  of  the 
coupling  constants  and  asymmetry  parameters  may  be  used  in 
conjunction  with  Eqs.  4.3-4  to  4.3-6  to  calculate  orbital 
populations  assuming  the  same  EFG  axes  as  in  HMX.  If  c^  is  held 
constant  at  l.lOe,  the  values  of  a  and  shown  in  Table  4.3-IV 
are  obtained.  aNlT  is  slightly  smaller  in  RDX  than  in  HMX 
indicating  that  the  amount  of  NN  ir-character  is  slightly  greater 
in  RDX.  This  conclusion  is  consistent  with  the  larger  value  found 
for  the  *^N  coupling  constant  in  RDX  compound  to  HMX.  ^4 
However,  the  value  of  b.,.,  is  found  to  be  less  in  RDX  than  HMX 

NN  - 

which  is  consistent  with  the  slightly  longer  average  N-N  bond  in 

o  AC  o 

RDX  (1.38  A)  compared  to  B-HMX  (1.36  A).  Therefore,  the  trend 
in  a„  compared  to  b..M  seems  to  be  contradictory  in  RDX,  yet  each 
value  finds  supiport  in  the  results  of  other  experiments.  It  is 
certainly  possible  that  the  limits  of  the  Townes-Dailey  model  have 
been  reached  as  far  as  the  NQR  data  are  concerned.  Additional 
work  on  this  problem  is  planned. 
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Molecular  Motion;  From  the  temperature  dependence  of  the 

coupling  constant  for  the  axial  and  equatorial  amine  nitrogen 

atoms,  one  can  deduce  information  about  the  nature  of  motion  in 

2 

HMX .  The  change  in  e  Qq/h  with  temperature  is  shown  in  Figure 

4.3-4.  The  often  found  slow  decrease  in  the  EFG  with  increasing 

temperature  occurs  in  HMX.  However,  the  temperature  coefficients, 
2 

d(e  Qq/h)/dT,  for  the  axial  and  equatorial  amine  sites  are 
different.  Three  temperature  regions  exist  where  the  coefficients 
at  each  site  are  reasonably  constant.  These  are  the  77-270° K, 
270-370° K,  and  370-420° K  regions.  Table  4.3-V  gives  the 
temperature  coefficient  values.  Figure  4.3-5  graphically 
illustrates  the  much  more  rapid  decrease  in  d(e  Qq/h)/dT  at  the 
equatorial  amine  sites  compared  with  the  axial  site  with  increase 
in  temperature.  No  fir  -  ord^r  phase  transitions  occur  in  HMX 
over  this  temperature  x^nge. 

2 

The  temperature  dependence  of  e  Qq/h  as  well  as  its 
difference  for  the  axial  and  equatorial  amines  originates  in  the 
temperature  dependence  of  the  librational  motions  of  the  molecule. 
These  motions  occur  about  the  X,Y  and  Z  EFG  axes  for  each  amine 
site.  When  the  torsional  motions  are  described  by  an  Einstein 
model,  the  root  mean  square  of  the  amplitudes  about  the  x,  y  and  z 
inertial  axes  of  HMX  are  expressed  as  a  function  of  temperature 
according  to  eq. (4.3-7). 


<e2>  = 

l 


47T  IiVi 


1/2  + 


exp(hv^/kT)  -1 


;  l  =  x,y, z 


(4.3-7) 
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TABLE  4.3-V 


THE  TEMPERATURE  COEFFICIENTS  FOR  THE  AXIAL  AND 
EQUATORIAL  AMINE  ATOMS  OF  HMX.  THE  PARENTHETICAL 
NUMBERS  ARE  THE  COEFFICIENTS  OF  CORRELATION  OVER 
THE  TEMPERATURE  REGION  GIVEN 


Temperature  Range.  °K 

d(e2Qq/h)/dT, 

— 

Hz/°K 

Equatorial 

Axial 

77-270 

-0 . 085 (-0 . 9S8) 

-0.1 27 (-0 . 993) 

270-370 

-0.181 (-0.990) 

-0.1 58( -0  985) 

370-426 

-0.365(-0.992) 

-0.262(-0.994) 
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Figure  4.3-5.  d(e  Qq/h)/dT  for  each  HMX  site  shows  greater  temperature 
sensitivity  at  the  Equatorial  site  than  at  the  Axial 


»  th 

is  the  effective  moment  of  inertia  about  the  i  axis  and 

is  the  librational  frequency.  To  the  extent  that  the  principal 

axes  of  the  inertial  moment  coincide  with  the  EFG  tensor  principal 

axes,  the  change  in  the  coupling  constant  and  asymmetry  parameter 

44  45 

with  temperature  can  be  expressed  in  Eqs.  4.3-8  and  4.3-9.  ' 


e2Qq/h  =  e2QqQ/h 


\  - 1 

f  2  2  1 

<0  >  +  <0  > 

+  1 

f  2  2  1 

<0^>  -  <0  > 

2 

x 

2 

y 

(4.3-8) 


n 


n  _e2Qq/h 
°  e2Qqo/h 


1  -  1/2 

f  2  2  1 

<0  >  +  <0  > 

+  f 

f  2  2  1 

<e^>  -  <eS 

-2<02> 

x  y 

.  j 

2n 

Y  ,  x 

l  ‘  A 

z 

(4.3-9) 


It  can  be  seen  from  Eq  4.3-8  that  librational  oscillations 

about  the  Z  axis  will  hardly  affect  the  coupling  constant. 

2 

Therefore,  e  Qq/h  is  mostly  a  function  of  librational  motion  about 

all  three  EFG  axes.  The  inertial  moment  axes  in  HMX  can  be 

calculated  by  using  the  atomic  coordinates  taken  from  neutron 

diffration  data  57  based  on  the  convention:  I  >1  ,  >1  .  The 

z  y  x 

values  of  I  in  HMX  are  found  to  be  I  =  149.7,  I  m  277.9,  and  I 

p  x  y  z 

=  374.2  in  units  of  amu  A-  .  The  origin  of  I  lies  at  the 
molecular  inversion  center.  The  EFG  axes  deduced  in  the 
preceding  section  do  not  precisely  coincide  with  these  inertial 
axes,  but  there  is  qualitative  coincidence  between  the  two  axis 
systems  shown  in  Figures  4.3-1  and  4.3-2. 


BLOCK 
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One  can  consider  the  z  inertial  axis  in  order  to  illustrate 

the  effect  on  the  EFG  of  motion  about  an  inertial  axes.  For  the 

axial  amines,  the  Z  axis  of  the  EFG  is  approximately  parallel  to 

the  z  inertial  axis,  while  the  X  and  Y  EFG  axes  are  approximately 

perpendicular  to  it.  For  the  equatorial  amines,  the  Y  EFG  axis 

approximately  coincides  with  the  z  inertial  axis  while  the  X  and  Z 

EFG  axes  are  perpendicular  to  it.  Therefore,  motion  about  the  z 

inertial  axis  causes  a  and  6  to  be  larger  than  0  for  the  axial 

x  y  z 

amines,  while  0  and  0  are  larger  than  0„  for  the  equatorial 
'  y  z  x 

amines.  9,0,  and  0  refer  here  to  librational  motion  about  the 
x  y  2 

X,  Y,  and  Z  EFG  axes.  By  using  Eqs.  (4.3-8)  and  (4.3-9)  one  can 

2 

reason  that  motion  about  the  z  inertial  axis  makes  d(e  Qq/h)/dT 


larger  for  the  axial  site  than  for  the  equatorial  site.  Similar 
reasoning  suggests  that  motion  about  the  y  inertial  axis  affects 
the  temperature  coefficient  of  the  axial  and  equatorial  sites  in 
approximately  the  same  way  while  motion  about  the  x  inertial  axis 
should  lead  to  [d (e2Qq/h) /dT] Equatorial  >  [d^Qq/hJ/dT]^.^. 


A  scheme  for  the  motion  in  HMX  as  the  temperature  increases 
toward  the  phase  transition  to  6-HMX  at  436°K  can  be  extracted 
from  these  data.  In  the  low  temperature  domain  (70-270°K)  the 
lattice  is  quite  rigid  but  some  librational  oscillation  about  the 
z  inertial  axis  is  consistent  with  the  observed  relative  values  of 
the  temperature  coefficients.  This  motion  does  not  represent 
rotation  of  the  molecule,  but  only  libration.  However,  it  is 
worthy  of  note  that  motion  about  this  axis  occurs  in  flat 
molecules  like  benezene&7  and  s-triazine68  which  bear  a 
resemblance  in  shape  to  0-HKX. 
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As  the  molecule  passes  through  the  room  temperature  region, 
the  importance  of  motion  about  the  x  (and  perhaps  y)  inertial  axes 
appears  to  become  increasingly  important  and  causes  the 
temperature  coefficient  of  the  equatorial  amines  gradually  to 
surpass  that  of  the  axial  amines.  In  the  higher  temperature 
domain  approaching  the  g-><S-HMX  phase  transition,  oscillations 
about  the  x  inertial  axis  appear  to  dominate. 

The  motion  about  both  the  x  and  y  axes  represents  oscillation 
about  the  pivot  formed  by  the  most  important  intermolecular 
electrostatic  interaction  in  crystalline  B-HMX:  the  contact 
between  an  oxygen  atom  of  a  neighboring  molecule  and  the  C^N2C2 
fragment  in  the  axial  portion  of  the  molecule  shown  in  Figure 
4.3-2.  This  interaction  leads  to  a  chain-like  structure  in  the 
6-HMX  lattice  (Figure  4.3-6)^  which  occurs  approximately  along 
the  z  inertial  axis  of  the  HMX  molecule.  Pivoting  of  the  molecule 
about  this  axis  interaction  would  yield  a  larger  temperature 
coefficient  in  the  low  temperature  region  for  the  axial  amines 
compared  to  the  equatorial  amines.  As  the  lattice  approaches  the 
phase  transition  to  <S -HMX,  libration  about  the  x  inertial  axis 
takes  place.  Like  the  z  axis  libration,  this  motion  involves  a 
type  of  pivoting  about  the  interaction  but  is  more  disruptive  of 
the  interaction.  Such  an  analysis  does  not  preclude  the 
possiblity  that  other  motions  in  HMX  are  important  as  well, 
particularly  because  the  ring  conformation  must  also  change  during 
the  phase  transition.  However,  disruption  of  the  most  important 
cohesive  interactions  is  a  prerequisite  to  the  B-+-6  transition  in 
HMX  and  this  is  found  to  take  place. 
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Figure  4.3-6.  A  view  along  the  "chains"  of  B-HMX  molecules 

which  form  as  a  result  of  the  interaction  shown 
in  Figure  4.3-3,  Oscillations  of  the  molecules 
along  this  chain  are  primarily  responsible  for 
the  changes  in  d(e2Qq/h)/dT  with  temperature. 
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The  results  presented  here  leave  little  doubt  that  the  x>st 
important  crystal  stablizing  interaction  in  0-HMX  both  in  terms  of 
electrostatic  energy63  and  molecular  dynamic  motion  is  the 
0...CNC  interaction  between  the  axial  N02  groups  in  neighboring 

HMX  molecules. 
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5.0  CONCLUSIONS  AND  RECOMMENDATIONS 


A  summary  of  the  specific  results  obtained  during  this 
program  is  given  in  Section  1.2.  Here  we  list  some  additional 
conclusions  we  reached  as  a  result  of  this  work: 


a)  Fourier  transform  and  pulse  NQR  techniques 
proved  to  be  sufficiently  sensitive  to  detect, 
with  adequate  resolution,  the  14N  NQR  spectrum 
of  TNT  at  all  temperatures  of  interest. 

b)  Efforts  to  detect  the  very  weak  TNT  lines  over 
the  last  five  years  resulted  in  the  development 
of  new  excitation  sequences,  i.e.  the  multiple 
pulse  SLSE  and  SORC  sequences,  and  in  improved 
FT  methods. 

o)  It  has  become  possible  to  detect  the  presence 
or  absence,  as  well  as  the  relative  abundance 
of  orthorhombic  and  monoclinic  TNT  in  a  given 
volume  of  space. 


Although  this  program  enjoyed  considerable  success, 
further  work  along  the  following  lines  would  still  be  useful 
to  a  more  nearly  complete  understanding  of  the  structure  and 
bonding  of  TNT: 


a)  Measurement  of  the  TNT  spectrum  at  liquid  helium 
temperatures  would  give  a  firmer  foundation  to 
the  determination  of  the  librational  motions  of 
the  TNT  nitro  groups.  Similarly,  to  study  the 
spectrum  and  the  relaxation  t^r.es  up  to  the 
melting  point  might  be  very  useful  in  understanding 
the  thermal  decomposition  of  the  molecule. 

b)  The  complex  (TNT)2*HN5  which  has  been  identified 
as  the  nucleating  agent  for  growth  of  the  stable 
(monoclinic)  phase  of  TNT  from  the  melt,  is 
currently  of  unknown  structure  and  conformation. 

An  NQR  study  of  structure  and  bonding  in  this 
system  would  be  invaluable,  and  would  complement 
any  eventual  x-ray  structure  determination. 


5-1 


BLOCK 


6 . 0  ACKNOWLEDGEMENTS 


The  authors  would  like  to  thank  our  technical  monitor 
Dr.  Charles  Boghosian,  for  his  support  during  the  tenure  of 
this  project  and  Dr.  Frank  Owens  for  his  continued  interest 
and  encouragement  during  the  course  of  this  work.  Dr.  Stanley 
M.  Klainer  is  acknowledged  as  the  man  without  whose  efforts 
the  original  proposal  would  never  have  been  written. 

Finally,  one  of  us  (RAM)  would  like  to  thank  the 
management  and  staff  of  Block  Engineering  for  the  warm 
hospitality  they  extended  to  him  while  he  served  as  principal 
investigator  for  this  project. 


6-1 


BLOCK 


APPENDIX  A 
NQR  THEORY 
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A. 1  INTRODUCTION 


The  technique  of  NQR  developed  as  an  outgrowth  of  the 
discovery  of  NMR  in  the  late  forties.  The  theory  of  nuclear 
quadrupole  interaction  in  solids  was  first  given  by  R.V.  Pound 
at  Harvard  in  1950,  and  later  that  year  Dehmelt  and  Kruger 

35 

in  Germany  observed  the  NQR  signal  from  the  chlorine-35  (Cl  ) 

nucleus  in  transdichloroethylene.  In  the  following  year 

Watkins  and  Pound  detected  the  first  quadruoole  resonance  from 
14 

Nitrogen-14  (  N)  in  HMT  (hexamethylenetetramine) ,  BrCN  and  ICN. 
Experimental  difficulties  delayed  the  utilization  of  the  NQR 

and  until  the  late  sixties  only  a  few  compounds  were  studied. 

14 

Two  chief  difficulties  arise  in  N: 

(1)  the  low  frequency  range  of  the  resonance  lines 
(0.5  to  5  MHz)  which  makes  them  weak  and  hard 
to  detect,  and 

(2)  the  unusually  strong  intermolecular  interactions 
which  are  characteristic  in  nitrogen-containing 
compounds  preclude  on  one  hand  the  proper  transfer 
of  radio-frequency  energy  to  the  crystalline  lattice, 
and  on  the  other  hand  make  the  electric  field 
gradient  very  sensitive  to  the  purity  and  crystal¬ 
linity  of  the  material. 

This  results  in  transitions  which  are  easily  saturated  and  lines 
which  are  too  broad  to  be  detected  with  conventional  (cw)  spec¬ 
trometers.  With  the  development  of  pulsed  methods  and  more 

sensitive  spectrometers  these  difficulties  can  be  overcome. 

14 

At  the  preset*-  time  N  NQR  is  successfully  being  applied 
using  the  Block  Engineering,  NQR/FFT  Spectrometer. 
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A. 2  BASIC  CONCEPTS  OF  NQR 


In  order  to  understand  the  origin  of  nuclear  quadrupole 
resonance,  we  may  visualize  a  nucleus  as  a  classical  distri¬ 
bution  of  positive  charges,  p.  (r)  ,  over  a  volume  of  character- 

N  -13 

istic  linear  dimensions  of  the  order  of  nuclear  radii:  10  cm 

(Figure  A-l) .  A  charge  cloud  extending  over  a  volume  of  linear 
dimensions  of  the  order  of  several  Angstroms  generates  an  electro¬ 
static  potential  <j>  (r)  which  can  be  considered  as  varying  slowly 
over  the  region  occupied  by  the  nucleus.  The  electrostatic 
energy  of  the  system  expressed  as 


W  = 


PN(r)<j>  (r)  dr 


(A-l) 


can  then  be  expanded  in  terms  of  the  moments  of  the  nuclear  charge 
distribution: 


W  =  c»>  ( 0 )  qN  +  p*2(0)  -  ~  l  +  ••• 

ij 


(A-2) 


where 


PN(r)d? 
PN(r)  r  dr 

pN(?)XiXjd? 


E 


-V<j) 

3E. /3x . 
i  3 


nuclear  charge 
dipole  moment  vector 
quadrupole  moment  vector 
electric  field  vector 
electric  field  gradient  tensor 
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From  quantum  mechanical  considerations  it  can  be  proven  that 
the  only  non-vanishing  terms  in  (2)  are  the  first  and  third  ones. 
Then,  if  the  energy  of  the  system  is  referred  to  the  constant 
first  term: 

WQ  =  “  [w  +<M0)qN 

0^^  is  diagonal  and  is  related  to  the  components  of  the  nuclear 
spin  1^  according  to 


=  ?  A 


(A- 3) 


eQ 

K2I-1) 


3Ii2  -  1(1+1) 


(A- 4) 


where  I 
moment . 


Since 


j 


h 


and  eQ  is  defined  as  the  nuclear  quadrupole 


^q^j  =  0  (Laplace  equation:  A<j>  =  0) 

the  quantum  mechanical  expression  of  the  quadrupole  energy 
(the  Hamiltonian)  takes  the  particular  simple  form: 


HQ  21(21-1)  |  ^ij*i  *  (A~5) 

which  can  also  be  expressed  in  the  following  way: 

KQ  -  [(3IZ2  -  l2»  +  "(Ix2  -  Iy2)]  ,A-6) 

2 

where  e  qQ  is  defined  as  the  quadrupole  coupling  constant  of 
the  system  and  n  is  defined  as  the  asymmetry  parameter  of  the 
electric  field  gradient  (EFG) : 
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where 


l<3zzi  >  !qyyi  >  lqxxl  by  convention- 

For  nitrogen-14,  of  nuclear  spin  1=1,  the  solution  of 
the  Schroedinger  equation  HQ  ipN  =  E^N,  where  ij>N  is  the  nuclear 
wave  function  can  be  shown  to  lead  to  a  three  level  system  of 
energies  given  by 

Ez  =  e2qQ/2 

Ex  =  V1_n)/2  (A-8) 

E  =  E  (1+n) /2 
y  £ 


Transitions  between  these  levels  can  be  induced  with  oscil¬ 
lating  magnetic  fields  of  the  proper  (resonant)  frequencies. 
The  frequencies  of  these  transitions  are  (see  Figure  A-2) 
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Conversely,  determining  any  two  of  the  three  NQR  frequencies 
completely  describes  the  EFG  in  the  vicinity  of  the  nitrogen 
nucleus: 


e2qQ 


j  (v+  +  v_) 


2<v+  -  v.) 


(A-10) 


2  _ 
e  qQ 


Experimental  techniques  for  the  determination  of  NQR 
spectra  can  be  divided  into  three  categories: 

(a)  Continuous-wave  (cw)  methods 

(b)  Transient  methods 

(c)  Double  resonance  techniques 

CW  methods  have  enjoyed  great  popularity  for  they  are 
simple  and  inexpensive.  The  limitations  of  these  techniques 
are  so  severe,  however,  that  at  least  in  the  Nitrogen-14 
NQR  field  they  are  rarely  used  any  more.  Transient  methods 
include  the  superregenera tive  technique,  which  ,  for  Nitrogen-14 
is  bad  and  rarely  used,  and  pulsed  methods  that  monitor  the 
resonance  of  one  type  of  nuclear  species  while  another  is 
being  perturbed.  Even  though  transient  methods  are  extremely  sensitive 
they  do  not  have  as  high  a  resolution  as  pulsed  FT  methods,  and 
are  not  the  .nethod  of  choice  for  this  particular  study. 

The  Nitrogen-14 ,  pulsed  techniques  that  operate  in  the 
so  called  spin-echo  mode  are  particularly  good  (for  other 
nuclei  the  reverse  holds  true,  and  CW  or  superregenerative 
methods  are  better) . 


BLOCK 


To  understand  the  reasons  for  this  advantage  let  us 
begin  by  defining  some  basic  magnitudes  that  describe  an 
NQR  line.  In  the  frequency  domain  an  NQR  line  is  fully 
described  by  its  frequency  and  by  a  normalized  shape 
function  S(w).  For  example: 


S(w) 


1  + 


w-wo 


Aw 


(A-ll) 


Spin-echo  techniques  operate  in  the  time  domain,  there 
are  three  parameters  that  fully  describe  an  NQR  absorption 
line 


(a)  Spin-lattice  relaxation  time(T^).  A  magnitude  which 
establishes  the  typical  "waiting"  time  to  observe  a 
resonance  a  second  time. 

(bj  Spin-spin  relaxation  time.  (T2 ) .  A  magnitude  wnich 

establishes  the  time  scale  for  the  observation  of  spin-echoes. 

(c)  A  spin-echo  shape  function  G(t).  The  "width"  of 

* 

this  function  is  defined  as  2T2.  G ( t )  is  the 

fourier  transform  of  the  line  shape  function 
*  1 

S (w) ,  thus  T2  a  ftw  * 

CW  and  superregenerative  methods  are  poor  techniques  for 
Nitrogen-14  because  of  the  generally  large  and  large  Aw 
values  associated  with  an  NQR  line  (to  avoid  saturation  in 
CW  the  level  has  to  be  reduced  to  very  small  values) .  Broad 
lines,  a  very  common  occurrence  with  Nitrogen-14,  are  weak 
since  their  areas  are  constant.  Spin-echo  methods  are,  on  the 
other  hand,  unaffected  by  long  T^  values  and  since  G(t)  and 
S(w)  are  related  by  a  Fourier  transformation  the  equality 
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S(w)  dw  =  1 


(A-12) 
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holds  (Figure  A-3) .  Thus,  broad  lines  in  the  frequency  domain  do 
not  affect  the  maximum  intensity  of  the  echo  signal.  In  addi¬ 
tion,  pulsed  methods  for  Nitrogen-14  are  particularly  advanta¬ 
geous  because  T2  times  are  generally  long  (of  the  order  of  mil¬ 
liseconds)  .  This  is  a  result  of  the  "spin-quenching"  effect 
for  n  ^  0  which  substantially  reduces  spin-spin  interactions. 


The  Nitrogen-14  NQR  spectrum  of  a  molecular  system  is 
determined  by  placing  about  25  grams  of  sample  inside  the 
inductor  of  a  resonant  circuit  which  is  then  subjected  to  a  series 
of  radio  frequency  pulses  of  frequency,  f.  Whenever  the  fre¬ 
quency  of  these  pulses  satisfies  the  resonance  condition 
f  =  v0  (where  vQ  is  one  of  the  quadrupole  frequencies)  absorp¬ 
tion  of  energy  takes  place  and  is  retransmitted  as  a  series 
of  signals  as  is  shown  schematically  in  Figure  A-4.  Therefore, 
by  monitoring  and  detecting  the  transmitted  signals  as  a  function 
of  the  frequency  of  the  pulses,  the  energy  levels  of  the 
quadrupole  nucleus  are  completely  determined. 


Since  a  narrow  RF  pulse  has  a  Fourier  spectrum  of 
finite  width,  the  pulse  frequencies  can  be  varied  in  discrete 
steps  (5-10  kH2)  so  as  to  cover  a  larger  frequency  range  for 
a  given  observation  time.  In  this  manner  the  total  frequency 
range  is  only  limited  by  the  "waiting  period"  between 
successive  observations.  This  waiting  period  depends  on  the 
spin  lattice  relaxation  time  (T^)  of  the  particular  compound 
which  is  being  investigated.  This  can  be  quite  long  (5-10  min.) 
for  Nitrogen-14  NQR. 


A- 9 


BLOCK 


FREQUENCY  DOMAIN 


TIME  DOMAIN 

FID 


Figure  A- 3 

Pulsed  KQR  Signals  and  Noaenclature 


Typical  N  Pulse  Sequence 
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B.  1  LIBRARY  OF  TNT  SPECTRA 


In  this  section  we  collect  a  library  of  TNT  spectra  to 
provide  a  record  of  the  Nitrogen-14  NQR  lineshapes.  The 
legend  of  each  plot  gives  the  crystal  prase  of  the  TNT  sample, 
the  temperature  (77°K)  and  the  number  of  coherently  added 
echo  signals  obtained  in  spin-locked  spin- echo  (SLSE) 
experiment . 

Each  plot  consists  of  the  time-domain  NQR  signal  displayed 
at  the  top,  and  the  modulus  Fourier  transform  of  the  echo 
signal  displayed  at  the  bottom.  Before  computation  of  the 
transform,  SORC-type  signals  seen  at  the  beginning  and  end  of 
the  time=domain  data  (see  Section  2.1)  were  apodized  by 
multiplication  with  a  "box  car"  function.  Thus,  only  the 
center  signal  (the  spin-echo)  was  transformed. 

The  caption  under  each  graph  lists: 

a)  The  date  on  which  the  plot  was  made 
and  a  plot  sequence  number 

b)  The  TNT  phase 

c)  The  reference  frequency  (zero  on  the 
plot  scale) 

d)  Sample  temperature 

e)  The  number  of  echos  coadded. 
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C.l  BROADBAND  FT  SPECTRUM  OF  EXCEPTIONALLY  WIDE  NQR  LINES 


In  Appendix  C  we  collect  the  broadband  FT  spectrum 
of  the  three  lines  in  the  hydrogen-bonded  complex  Hexamethy¬ 
lenetetramine  Hexahydrate,  HMT’EH^O.  The  experiment  methods 
are  described  in  Section  2.2.1.  We  shall  limit  the  discussion 
in  this  Appendix  to  the  actual  data  obtained;  an  analysis 
of  the  lineshape  is  beyond  the  scope  of  this  work. 

Figures  C-i,  C-2,  and  C-3  are  the  line-shapes  of  the 
three  NQR  lines  in  HMT*6H20,  respectively  called  the  v+,v_ 
and  vq  lines.  These  spectra  were  obtained  by  adding  together 
two  "normalized"  spectra  obtained  by  setting  the  spectrometer 
frequency  above  and  then  below  the  average  resonance  frequency. 
The  normalization  was  done  by  dividing  the  raw  spectrum  by 
the  frequency  response  of  the  sample  coil  and  receiver,  and 
.is  not  explicitly  shown  here.  Figures  C-4,  C-5,  and  C-6 
are  composites  showing  the  raw  and  normalized  spectra  from 
which  Figures  C-l,  C-2  and  C-3  were  obtained. 
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Figure  C-l.  Line  Shape  for  the  v+  line  in 

HMT  ‘6H2O  at  77°K,  The  abscissa 
shows  increasing  frequency  from 
3235.0  kHz. 
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Figure  C-2.  Line  Shape  for  the  v_  Line  in 

HMT'SH^O  at  77 °K.  The  abscissa 
shows  increasing  frequency  from 
3140.0  kHz. 


BLOCK 

t*)C*+44Ca***Q 


1(0  20  30  4©  50 

FREQUENCY  CKHZ7 

NORMALIZED  HMT*SH2Q  ZERO  LINE  CA+F8>  RATE  l.S  77K 

Figure  C-3.  Line  Shape  for  the  v  line  in 

HMT^GH^O  at  77°K.  The  abscissa 
shows  increasing  froquency  from 
3560.0  kHz. 
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Figure  C-4 


a)  Time-domain  spin-echo  signal  and  modulus  Fourier 
transform  for  the  v+line  of  HMT’S^O  at  77°K. 

The  abscissa  shows  frequency  increasing  from  the 
transmitter  frequency  at  3225  kHz. 

b)  Spectrum  in  (a)  corrected  for  sample  coil  Q  and 
receiver  response. 

c)  Same  as  (a)  but  abscissa  shows  frequency  decreasing 
from  the  transmitter  frequency  at  3255  kHz. 

d)  Spectrum  in  (c)  corrected  for  sample  coil  Q  and 
receiver  response. 
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Figure  C-5 


a)  Time-domain  spin-echo  signal  and  modulus  Fourier 
transform  for  the  v_  line  of  HMT*6H20  at  77°K. 

The  abscissa  shows  frequency  increasing  to  the 
right  from  the  transmitter  frequency  at  3560  kHz. 

b)  Spectrum  in  (a)  corrected  for  sample  coil  Q  and 
receiver  response. 

c)  Same  as  (a)  but  abscissa  shows  frequency  decreasing 
from  the  transmitter  frequency  at  3590  kHz. 

d)  Spectrum  in  (c)  corrected  for  the  sample  coil  Q 
and  receiver  response. 
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Figure  C-6 


a)  Time-domain  spin-echo  signal  and  modulus  Fourier 
transform  for  the v0  line  of  HMT*6H20  at  77°K.  The 
abscissa  shows  frequency  increasing  from  the 
transmitter  frequency  at  3560  kHz. 

b)  Spectrum  in  (a)  corrected  for  sample  coil  Q 
and  receiver  response. 

c)  Same  as  (a)  but  abscissa  shows  frequency  decreasing 
from  the  transmitter  frequency  at  3590  kHz. 

d)  Spectrum  in  (c)  corrected  for  sample  coil  Q  and 
receiver  response. 
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PRINCETON  APPLIED  RESEARCH  INTERFACE 
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In  the  original  Block  NQR  spectrometer  system,  the  analog 
signal  with  the  associated  status  and  clock  timing  information  was 
bussed  through  a  screen  room  wall  via  optocouplers.  The 
signal  was  subsequently  digitized  and  coadded  by  a  Digilab® 

FTS  14  data  system  in  real  time.^ 

The  maximum  clock  rate  with  this  arragement  was 
50  kHz,  which  limited  the  bandwidth  to  25  kHz. 

In  order  to  extend  this,  a  Princeton  Applied  Research  (PAR) 
signal  averager  Model  4202  was  used  to  do  the  coaddition. 

This  had  clock  rates  to  200  kHz,  and  was  directly  usable 
on  the  NQR  instrument  with  the  available  output  gates. 

However,  it  was  necessary  to  read  out  the  data  from 
the  averager  into  the  data  system  for  further  processing 
in  some  manner. 

The  method  used  was  to  build  an  interface  which 
read  out  the  digitized  data  from  the  averager,  converted 
it  to  analog  form,  and  added  a  clock  and  status  output. 

This  data  was  then  bussed  through  the  optocouplers  and 
coadded  in  the  data  system.  The  advantage  of  this 
approach  was  that  it  was  easily  imp] ementable  and  mini¬ 
mized  hardware  changes.  Also,  the  PAR  4202  has  provisions 

2 

for  a  hook  up  in  this  manner. 

*A  Pulsed  NQR-FFT  Spectrometer  for  Nitrogen-14,  J.C.  Harding,  Jr. 
et  al.  Journal  of  Mag.  Resonance  36,  21-33  (1979). 

^ ‘Model  4202  Signal  Averager  Operating  and  Service  Manual 
Princeton  Applied  Research  1978. 
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T>e  circuit  shown  in  Figure  D-l  performs  these 
functions.  Control  of  the  PAR  4202  digital  output  is 
performed  by  the  signals  DIG  ORIGIN  RESTORE,  DIG  DATA  STROBE 

- 3 

and  DIG  DATA  VALID.  The  first  signal  resets  the  origin, 
while  the  other  two  signals  handshake  successive  points  out. 

The  circuitry  is  fairly  straight  forward.  Switching 
from  stop  to  start  initializes  the  logic  and  the  PAR  4202. 

Subsequently  the  status  is  enabled  and  the  74LS629 
clock  starts.  For  each  clock  point  the  DIG  OUTPUT  STROBE, 
DIG  DATA  VALID  signals  cause  the  PAR  4202  to  put  the  digital 
value  of  the  data  point  on  its  digital  output  lines.  The 
monostable  at  U7  allows  the  D/A  converter  to  settle,  and 
then  enables  the  clock  out  to  the  optocouplers. 

At  the  2048th  clock  point,  if  the  start  switch  is 
still  in  the  start  position,  the  logic  is  reinitialized, 
and  the  sequence  begins  again;  otherwise  it  stops. 

Each  data  word  of  the  PAR  4202  is  28  parallel  bits, 
while  the  digital  to  analog  converter  is  only  12  bits. 

In  order  to  properly  choose  the  appropriate  12  bits  for  the 
converter  in  a  reasonable  manner,  the  28  bits  were  bussed 
to  28  successive  terminals  and  a  connector  consisting  of 
12  successive  matching  terminals  was  constructed.  The 
choice  of  which  successive  12  bits  to  use  is  done  by 
manually  moving  the  connector  and  monitoring  the  output 
on  a  scope. 

^Model  4202  Signal  Averager  Operating  and  Service  Manual 
Princeton  Applied  Research  1978. 
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